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1. Introduction 

It is now widely believed that the 4 -scalar theory in four space-time dimensions is 
trivial. Consider the Lagrangian density 

( j> 4 (1) 

relevant to the scalar sector of standard model (SM). The evolution of the renormalized 
coupling constant A(M) at the scale M for large M, in the one loop approximation [1] is 
given by the equation 

(2) 



A(M) A(m) 167T 2 \m 

where m is some low energy reference scale. From (2), it is clear that A(M) diverges at the 
scale M = AL, where 

A L -mexp(167r 2 /3A(m)). (3) 

This is the 'Landau catastrophe'. The occurence of Landau catastrophe at any finite 
energy can thus be avoided only if A(m) = 0, which is the standard triviality scenario. 
One can however regard the Lagrangian (1) as an effective one; for a given A(m), the 
Lagrangian (1) then is a good description of the theory at most up to the energy scale A L 
as given by (3). The coupling constant A(m) can be related to Higgs mass M H . It is 
reasonable to expect that for a consistent effective theory M H should not exceed the scale 
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Abstract. A new exact closed form solution of Einstein's field equations is reported describing 
the space-time in the interior of a fluid sphere in equilibrium. The physical 3-space, t = constant of 
its space-time has the geometry of a 3-pseudo spheroid. The suitability of this solution for 
describing the model of a relativistic superdense star is discussed and the stability of the model 
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1. Introduction 

The non-linear nature of the Einstein's field equations is a consequence of the self- 
interaction of the gravitational field. This makes it difficult to obtain relativistic models of 
spherical stars based on exact solution of Einstein's field equations describing spherical 
distributions of matter. The actual properties in the central region of a relativistic compact 
star are not precisely known and so assumptions of general nature to obtain exact 
solutions of Einstein's field equations become necessary. It is desired that solutions 
should be physically plausible and at the same time simple in form. This problem has 
been considered by several authors [1-6]. 

Here we have investigated the gravitational significance of space-times whose physical 
space obtained as t = constant section has the geometry of 3-pseudo spheroid and it 
is shown that it can describe spherical compact distributions of matter. The form of 
the space-time metric and its general features are discussed in 2. It is shown that the 
geometry of the space-time is governed by two parameters R and K. The space-time 
metric corresponding to a new exact closed form solution of Einstein's field equations is 
written in 3. The suitability of this solution for describing the model of a superdense star 
is explored following the approach of [4, 5, 9] in the subsequent sections. 



2. Static pseudo spheroidal space-time 



dff 2 = dx 2 + dy 2 + dz 2 + dw 2 , (1) 

will have the cartesian equation 

w 2 x 2 +y 2 + z 2 _ , . 

.P" ~R 2 ~ ' () 

where b and R are constants. The section w = constant of the 3-pseudo spheroid are 
pseudo spheres, while sections jc = constant, y = constant and z = constant represent, 
respectively hyperboloids of two sheets. 
The parametrization 

Jt = rsin0cos<?!>, y = rsin0sin0, z=rcos0, w = b(l+ r 2 /R 2 } l/2 (3) 
of the 3-pseudo spheroid leads to 



da 2 = ~~ dr 2 + r 2 (d9 2 + sin 2 9d^ (4) 

where 

K=l+b 2 /R 2 . (5) 

The pseudo spheroidal 3-space given by eq. (4) is spherically symmetric and regular 
for K > I. It is flat when K = 1 and generates into open hyperboloid when K = 0. 
Following the Vaidya-Tikekar [4] approach we consider the space-time with metric 

d/ = e"Wdf 2 - 



Various aspects of spherical distributions of perfect fluid in equilibrium described by 
space-times with metric in a form similar to (6) have been investigated by Buchdahl [1]. 
Vaidya and Tikekar [4] have examined the geometrical features of the 3 -dimensional 
physical spaces obtained as / = constant hypersurfaces of such space-times. The space- 
time of the specific solution discussed by Buchdahl has the geometry of a 3-spheroid 
immersed in a 4-dimensional Euclidean space of Vaidya-Tikekar type. An extensive 
study of such solutions has also been done by Maharaj and Leach [6]. The 3-space of the 
metric (6) is characterized by pseudo spheroidal geometry and therefore the specific class 
of solutions in this set up is clearly distinct from the specific class of solutions reported by 
Buchdahl, Vaidya and Tikekar [4] and Maharaj and Leach [6]. 

3. Matter distribution on pseudo spheroidal space-time 

In this section the gravitational significance of a static, spherical distribution of matter in 
the form of perfect fluid is explored on the background of the space-time of the metirc (6) 
with energy-momentum tensor 

Ttj = ( P + p/c 2 }uiUj - (p/c 2 } gij . (7) 

Here p,p and u l respectively denote the matter density, fluid pressure and the unit 4- velocity 
field of the fluid. Since the fluid distribution is at rest, 

ii / = (0,0,0,e-"/ 2 ). (8) 

Q6 Pramana - T Phvs . Vnl $(\ TVn. 7 Fphrnarv 1QQ8 



Einstein's field equations 

_, I, 8?rG . 

Kij--Kgij = -z-Tij (9) 

for a specific choice of the curvature parameter K = 2 leads to the space-time metric 
written explicitly as 



ds 2 = AX/! + r 2 /R 2 + Bx(r] --^/\+ 2r 2 /R 2 dt 2 

V2 / 



sin 2 (W0 2 ), (10) 

where 

X(r] = -/I + r 2 //? 2 ln[\/2 v'l + >" 2 /# 2 + >/l + 2r 2 //? 2 ] (11) 



and A, B are arbitrary constants. 
The matter density and fluid pressure for the distribution (10) are expressed as 

8?rG 3 / 2r 2 \ / 2r 



(13) 

It is evident from eq. (12) that density is positive throughout the distribution. The 
gradient of density dp/dr is found to be negative indicating that p is decreasing radially 
outward. 



4. Size of the fluid sphere 

The total mass and size of the configuration can be estimated using the scheme given by 
Vaidya and Tikekar [4], as follows. 

Equation (12) determine p at the boundary r = a of the distribution as 

3(1+2^/3^) 

( 



We introduce the density variation parameter 

p(g)_ l+2a 2 /3R 2 

/\ 7f\\ "~"~ "~*) ) \ / 

where p(0) is the density at the centre. Since p is a decreasing function of r, A < 1 . 
Solving (15) as quadratic in (a 2 /R 2 } one finds 



(16) 



R 2 12A 
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The algebraic root assigning negative values to a 2 //? 2 is rejected to ensure that a/R is 
real. 

Equation (12) implies that the matter density at the centre is explicitly related with 
curvature parameter R as 

<W-- 07) 



Equation (17) therefore determines R in terms of p(a) and A. Equation (16) then 
determines the boundary radius a of the distribution. Thus the size of the configuration is 
determined in terms of the surface density p(a] and density variation parameter A. 

5. Physical plausibility 

Any physically acceptable solution must comply with the following conditions: 

(i) The matter density p and fluid pressure p should be non-negative throughout the 

distribution. 

(ii) The gradients dp/dr and dp/dr should be negative, 
(iii) The speed of sound should not exceed the speed of light as implication of causality 

fulfillment. 

(iv) The interior metric should match continuously with the Schwarzschild exterior 
solution 



(18) 



at the boundary surface r a of the distribution where p(d] 0. 
The continuity of the metric coefficients give 



and 



V/l - 2m/a = A-v/1 + a 2 /R 2 + B (X(a) - (l/N/2) v/1 + 2fl 2 //? 2 ), (20) 
where 



X(a] = v 'l+fl a /* 2 Mv2v'l + * 2 /* 2 + v / l + 20 2 /* 2 )- (21) 

The continuity of pressure across r = a requires that pressure to vanish on the 
boundary implying that 



= -B(X(a) + (1/2)0+ 2a 2 /R 2 }. (22) 

The constants A and B are determined from eqs (19), (20) and (22) as 



, 

' ( } 



B = _ . (24) 

\/2(H-2a 2 //? 2 )' 

while the total mass m is determined by (19). 

The expressions (23) and (24) for A and B when substituted in (13), one can find after a 
lengthy but straightforward computation that p > throughout the sphere. Further using 
the TOY equation we found that the pressure is decreasing radially outward. 

Using the expressions (12) and (13) for density and pressure and the values of A and B 
in the TOY equation it can be readily seen that the speed of sound will not exceed the 
speed of light in the central region and boundary of the distribution, ensuring the 
fulfilment of causality requirements. 

The expression for dp/dp, which represents the speed of sound in isentropic fluids, 
takes the form 



dp = irGR 2 (p + p/c 2 }(\ + 2r 2 //? 2 ) 3 [1 + (%-nGR 2 p/c 4 }(\ + 2r 2 /R 2 }} 

dp~ . (5 + 2r 2 //? 2 )(l+r 2 //? 2 ) ' ( ' 

At the centre, dp /dp has the value 

-dp\ = 7rG/? 2 [^(0) +p(0)/c 2 ](l + S7rGp(Q)R 2 /c*} 

Vo 5 

From eqs (12), (13), (23) and (24) we can show that 

p(0) - 3/?(0)/c 2 > (27) 

at the centre and using (17) and (27) it readily follows from (26) that 



0.2c 2 . (28) 



At the boundary, using (14), we have the expression 
3(1 + 2 2 /* 2 )(l + 2a 2 /3R 2 ) i 



dp) s 8(5 + 2a 2 /R 2 )(l + a 2 /R 2 } 



(29) 



The variation of dp /dp, which represents the speed of sound in isentropic fluids, is 
examined using numerical procedures for certain specific models in this set up. It is found 
that p/p is decreasing radially outward indicating that d(p/p)/dp > 0. However the 
speed of sound is found to be increasing radially outward for a number of models, which 
is an unsatisfactory feature of these solutions in view of the expectation that it should be 
decreasing radially outward which follows from equations of state found in literature. 
Since definite information about the equation of state for matter in nuclear density ranges 
is lacking, as argued by Knutsen [9], one must be careful in this respect. 

6. Superdense star model 

When all thermonuclear sources of energy are exhausted, a star will gradually cool down 
and in this process it will collapse gravitationally and form a compact star - white dwarf, 
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Table 1. Masses and equilibrium radii of superdense star models corresponding to 

K = 2 and p(d) = 2 x 10 14 gm/cm 3 . 



A 


tf(km) 


a (km) 


m(km) 


m/Mo 


A 


B 


0.900 


26.928 


04.875 


0.075 


0.050 


1.094 


-0.674 


0.800 


25.388 


06.828 


0.215 


0.146 


1.062 


-0.639 


0.700 


23.748 


08.271 


0.403 


0.273 


1.027 


-0.602 


0.600 


21.987 


09.429 


0.634 


0.429 


0.987 


-0.562 


0.500 


20.071 


10.382 


0.904 


0.613 


0.943 


-0.518 


0.400 


17.952 


11.162 


1.216 


0.825 


0.890 


-0.469 


0.300 


15.547 


11.771 


1.571 


1.065 


0.826 


-0.413 


0.200 


12.694 


12.176 


1.972 


1.337 


0.742 


-0.345 


0.125 


10.035 


12.291 


2.304 


1.562 


0.653 


-0.279 


0.100 


08.976 


12.274 


2.421 


1.641 


0.613 


-0.252 



Note: 1 Mo =1.475 km. 

neutron star, or black hole. Generally a star keeps its equilibrium with outward pressure 
against the self gravitational force. The model that we have presented here can describe 
the hydrostatic equilibrium conditions in such a superdense star with densities in the 
range 10 14 -10 16 gm/cm 3 . We take the matter density on the boundary r = a of the star as 
p(a) = 2 x 10 14 gm/cm 3 . Choosing different values for A, we determine the boundary 
radius a of the star and its total mass m in accordance with the scheme of 4 and 5. The 
value of m obtained is in kilometers. The mass of the star in grams is obtained using 
M = mc 2 /G. The results of these computations together with the values of constants A 
and B as determined by eqs (23) and (24) are given in table 1 . 

For A < 0.7 in table 1 , we get a set of physically viable models of relativistic compact 
stars. The models with A > 0.7 in the table have their equilibrium radii much smaller than 
that of a neutron star. The equilibrium models presented here take lesser values for radii a 
and total mass m, than the corresponding values given by Vaidya and Tikekar [4], 



7. Dynamic stability 

A sufficient condition for the dynamic stability of a spherically symmetric distribution of 
matter under small radial adiabatic perturbations has been developed by Chandrasekhar 
[7]. A normal mode of radial oscillation for an equilibrium configuration i.e. 

6r = (r) exp(iwr) (30) 

is stable when the frequency of oscillation u> is real and unstable when uj is imaginary. 
Chandrasekhar's [7] pulsation equation for the line element (6) is given by 



a; 2 



boundary 

centre 
boundary 



exp 



( + 



dr 



r dr 



dp fdu\ 
~ 



>dr, 



(31) 



where 

1 i O.-2 /2 

with K = 2. (32) 



The boundary condition to be satisfied at r = a is that the Lagrangian change in 
pressure 



where 7 is the adiabatic index. We must have 

~ = at r = a. (33) 

dr 

Following the method of Bardeen et al [8] and used by Knutsen [9] to investigate the 
stability of Vaidya-Tikekar models, we choose 



as the trial function where the new variable x is taken as x = 2r 2 //? 2 . 
The boundary condition du/dr = at r = a implies 



where b = 2a 2 /R 2 . 

The pulsation equation (31) for the metric (6) now takes the form 



+ \} 2 [L(x) +p(x) + q(x}} ' 
-8[L(*) +/>(*)] 



q(x}] ' 

)(* + 2) 



(34) 



(35) 



/-boundary /Sa + iA ( n+D}u 2 f b 

J / exp( ^p Y+P2L dr = / Rl 

Jcentre \ ^ / r JQ 

(36) 
where 

x + 2 l(x] - 
+ q(x] 



Table 2. The values of the integral on the right hand side of the 
pulsation equation (36) for some specific choices of the constants 
aj and b\ and A = 0.4. 



a\ 


bi 




Integral 


0.000 


-0.717 




1.476 


-0.776 


0.000 




0.536 


-0.858 


1.000 




0.084 


1.000 


-1.641 




3.447 


5.000 x 10 2 


-4.627 x 


10 2 


1.076 x 10 5 


-5.410 x 10 2 


5.000 x 


10 2 


1.246 x 10 5 


1.000 x 10 5 


-9.240 x 


10 4 


4.272 x 10 9 


-l.OOOx 10 s 


9.240 x 


10 4 


4.272 x 10 9 




B^=- 



We have evaluated the integral on the right of the equation (36) numerically for 
different values of b. It is found that the integral admits positive values for 0.243 < b < 
1.146 (i.e. 0.3 < A < 0.7) and for different choices small, large, postive or negative 
values of the constants a\,b\. 

We have reported in table 2 the values of the integral on the right side of the equation 
(36) evaluated numerically for certain specific choices of a\ and b\ for the model with 
A = 0.4 of table 1. 

This analysis indicates that these models with 0.3 < A < 0.7 will be stable. The space- 
time with pseudo spheroidal geometry for its spatial sections t = constant thus may admit 
the possibilities of describing interiors of superdense fluid stars in equilibrium. 
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1. Introduction 

Introduction of viscosity into the material content of model universes is gaining 
importance [1]. This is in the hope of having realistic description of the real universe 
mainly during the early stages of its evolution. Several workers have considered viscous 
fluid but only a few of them have tried to make it thermodynamically consistent. 
Thermodynamical consistency of the viscous fluid is undoubtedly essential as it 
leads to more realistic description. Moreover, if this consistency occurs in a way 
which allows the dissipation of viscosity causally, the reality of description is further 
enhanced. 

Initial efforts to obtain the above consistency were based on inadequate theories given 
by Eckart [2], and Landau and Lifshitz [3] - inadequate in the sense that these (theories) 
allowed viscous dissipation with super luminal speeds [4, 5]. However, this drawback was 
later on overcome by resorting to extended irreversible thermodynamics (KIT) which 
provided the relativistic causality in a natural way [6, 7]. Thus the viscosity obeying HIT 
may be a reasonable resort. 

In the frame work of standard (Robertson-Walker) cosmology, models have been 
derived which are causally well behaved as regards their viscous-fluid contents. Hiscock 
and Salmonsan [8] have considered the flat case with bulk viscosity. Zakari and Jou [9] 
and Maartens [10], while discussing it, have investigated the possibility of exponential 
inflation. Banerjee and Beesham [11] have obtained a power law solution as well. In the 

ri r\ _-f! o f roo \ITCH Vioirja r\-P /-r\nfO\ +Via o/-\lnti r\r\ rri\rt*in \\\r \ftir\At* 1 ? *^-nA T3at7V\f\ PI 01 V\iif -in +V\*a 



We consider the Robertson-Walker space-time given by 

ds 2 = -dt 2 + R 2 (t)[(l - kr 2 r'dr 2 + r 2 (d0 2 + sin 2 0dy? 2 )], (1) 

where k = 0, 1 and the energy momentum tensor Ty in the Einstein's field equation 
dj = Ty (for 87rG = 1 = c) (2) 

of the form 

T V = (P + Peff )// + Pe 



where p e ff = /? + II; p, p, II and U[ being respectively the energy density, equilibrium 
pressure, bulk viscous stress and the unit flow vector of the fluid co-moving orthogonally 
to the surface of homogeneity. We also consider the laws of 

1. Energy conservation 



2. Number conservation 

N l . = 0, N* = nu*. (5) 

alongwith 

3. The H-theorem 

^>0. (6) 

4. The Gibb's equation 



Tds = d(p/n)+pd(l/n), (7) 

where 5, n, T > 0, > and r > are respectively the specific entropy, number 
density, temperature, bulk viscosity coefficient and the relaxation coefficient for bulk 
viscous effects. A semicolon in the above indicates covariant differentiation. 

Equations (3)-(7) give rise to evolution equation [10, 11] for the bulk viscosity 

H + rli = -3H - (e/2)rn(3# + f/r - / - t/r), (8) 

where 



9 being the expansion scalar and an overhead dot standing for differentiation with respect 
to time t. Equations (l)-(3) lead to 

(R/R) 2 +k/R 2 = p/3 (9) 

and 
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Thus we get only three independent equations (8)-(10) in seven unknowns p, p, II, , r, T 
and R. So, as usual, we assume the following ad-hoc relations 

p = (7-l)p, 1<7<2, (11) 

= ap, a>0, <?>0, (12) 

r = t/p, (13) 

and 

r = /3p r , /3>0, r>0. (14) 

Now with the help of (9), (10) and (11), eq. (8) can be written in the form 
rH + (3/2)r(e + 2j)HH + H+ (l/2)er(f/r - / - T/T)H 
+ (9/4)e 7 r# 3 + (3/4)e 7 r(f/r - / - t/T 1 )/* 2 + (3/2) 7 // 2 



x [1/T + (3/2)eff - 2/0? + (c/2)(f/r - fa - T /T)] = 0. (15) 

Using (12)-(14), eq. (15) becomes 

H + (3/2) (27 + e}HH + a~ l [3(H 2 
x[H + 3~/H 2 /2 + (k/2R 2 )(3j - 2)] 
+ p/2# 2 )(3 7 - 2)(3e/2 - 2) - (9/2) (H 2 + kjR 2 }]H 
- (r +l)eH-(H- k/R 2 }/(H 2 + k/R 2 } 
x[H + (3/2)jH 2 + (*/2# 2 )(3 7 - 2)] = 0. (16) 



A general solution of eq. (16) is difficult to obtain. However, for the case k = 0, Maartens 
[10] has obtained an exponential solution whereas Banerjee and Beesham [11] have 
derived a power law solution. We notice in the present paper that R = t with q = 1/2 is 
also a solution provided: 

(i) Fork= 1: 

[V6(3 7 - 2)1(5 + 67) > when e = 

& = < 2V6(3 7 - 2)/4(5 + r) - 3 7 (1 + 2r) > 0, 

[ 0<r< (20-3 7 )/2(3 7 -2) whene=l. 

(ii) For k = 0: 

(V3(3j - 2)1(5 + 67) > when e = 



3. Discussion 



0, 
0<r< (20-3 7 )/2(3 7 -2) whene=l. 



For k = 1, the solution R = t, holds identically, though eq. (16) is not the proper way to 
go through in this case. Actually, in this case, p = = p. Thus no matter exists and the 
universe is empty. Thus R t is inconsistent when k = 1. 
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When k = 1, we obtain for the solutions corresponding to both the truncated (e = 0) 
and full (e = 1) theories the following expressions for concerned quantities: 

H = 2(2 -3 7 ) A 2 , 



Clearly, for / > 0, each of the quantities p, p, , r and T is positive, while II is negative. 
Moreover, both p and p tend to oo as t 0, and they tend to as t oo. Also 
= 3/t > oo and respectively when f > and oo. Thus the model starts with a big 
bang at ? = and its rate of expansion vanishes asymptotically. The temperature 
decreases gradually to zero when the expansion stops. We observe that II = p + p/3 
which restricts the model to be inflationary (actually, the condition for inflation R > 
requires, by virtue of (10), that -p e ff > P/3). The behaviour of the models are similar to 
that corresponding to the solutions for k 0, in which case the expressions for various 
quantities differ slightly by some positive constants only. The solutions for k 0, form 
the particular cases of that obtained by Banerjee and Beesham [11]. 
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Abstract. Particle production in curved spacetime has been discussed through the method of 
complex time WKB approximation. We consider Dirac equation in non-flat spacetime to 
understand particle production as particle-antiparticle rotation. The method is also generalized to 
understand particle production through parametric resonance. To understand the method of CWKB 
we consider particle production in Kasner spacetime as an example. 
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1. Introduction 

The Dirac equation in curved spacetime is of considerable interest in astrophysics and 
cosmology. Previous investigations in this direction are due to Fock [1], Tetrode [2], 
Schrodinger [3] and McVittie [4]. A general discussion of the m = neutrino case was 
given by Brill and Wheeler [5]. 

In astrophysics and cosmology the pair creation in expanding universes has been of 
considerable interest for several reasons. The choice of quantum vacuum is a very 
delicate object in curved spacetime and it greatly affects the calculations leading to pair 
creation. The works in this direction are due Chimento and Mollerach [6], Audretsch and 
Schafer [7], Castagnino and Mazzitelli [8], Biswas and Guha [9], Biswas et al [10], Guha 
etal[ll] and Biswas etal[l2]. Calculation of pair creation and others related to has also 
been carried by Lotze [13], Barut and Duru [14], Biswas et al [15], Sahni [16], Barut and 
Singh [17] and many .others. 

In astrophysics the study of Dirac equation has gained importance in strange matter 
calculations as well as in neutron star calculations [18], with and without magnetic field. In 
cosmology the effect of particle production has some physical considerations, mainly in 
the study of damping and avoidance of initial singularity. A study of pair creation of spin ^ 
particles (massive and massless) in Robertson-Walker universes was made by Parker [19], 
Lotze [13] and recently by Biswas et al [11, 12]. In the so called back-reaction problem, 
the effect of particle production back on the metric has not been incorporated successfully 
so far. Fischetti, Hartle and Hu [20], Anderson [21] tried to solve the back-reaction 
problem but the effect of particle production is not directly involved in the calculations. 



problem, but the calculations, though interesting, lack to project the dynamical evolution 
of spacetime that can only be solved through the solution of Einstein field equations with a 
term for particle production density in the right-hand side of Einstein equations. 

Various authors have studied particle production both in flat or non-flat spacetime. The 
common method is to define in and out vacuum and find the Bogolubov coefficients. The 
temporal equation resulting from Dirac equation in a Robertson-Walker framework has 
been found to differ [16, 17]. It has therefore been necessary to restudy Dirac equation in 
curved spacetime and find the correct temporal equation. It should be pointed out that the 
form of temporal equation is crucial in determining the particle production in expanding 
spacetime. In our previous works we followed a different procedure to calculate the pair 
production amplitude both in flat or non-flat R-W spacetime. In our approach [9-12], the 
turning points determined from the second order temporal Dirac equation are crucial in 
calculating the pair production amplitude. The turning points are fixed by the temporal 
equation in a given spacetime. To find out the correct turning points we investigate the 
derivation again without taking any results from other works. In this work we find that the 
Dirac equation in diagonal tetrad differs significantly from that of Barut and Singh [17]. We 
carry out the explicit derivation and find that our calculation is consistent with Sahni [16]. 

The present work, though a recheck of some previous calculations, also deals with 
some new results e.g., use of two dimensional temporal Dirac equation to understand 
particle production as particle-antiparticle rotation and parametric resonance particle 
production as a generalization of our result. 

In 2 and 3 we consider Dirac equation in expanding non-flat spacetime and reduce it 
to two dimensional form to study particle-antiparticle rotation. Section 4 deals with 
parametric resonance particle production and 5 deals with particle production in Kasner 
spacetime. We end up with discussion in 6. 

Recently particle production has been an active area of research [23-26] in cosmology 
to understand reheating mechanism after inflation, and the origin of cold dark matter. The 
method used [23-26] is parametric resonance. Our approach is an alternative one very 
similar to parametric resonance particle production. It is the aim of this work to highlight 
our approach with some specific examples. Though much work has been done on particle 
production, the present attempt sees the subject from a different approach with 
applicability in other directions. 

2. The Dirac equation in expanding non-fiat spacetime 

The Dirac equation in curved spacetime is 

= m^ (1) 



Here ^(x) are curvature dependent Dirac matrices, T^x] the spin connections and m the 
mass of the particle. 

2.1 The metric 

We start with Robertson-Walker type metric that describes a homogeneous isotropic non- 
flat universe given by 
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Particle production in curved spacetime 



ds 2 = dr 2 - a 2 (t) 



where k = +1, 0, 1, and for fc = this metric reduces to the conformally flat case 
order to understand the difference with flat-space equation we convert eq. (2) using 
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2.2 Affine connections 

Using g MJ/ and g^ 1 ' as given in equations (5) and (6) we calculate the affine connecti 
T with a, fj,, v = 0, 1 , 2, 3 to take the form 
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2.3 Curvature dependent Dirac matrices 

Henceforth curvature dependent Dirac matrices are written as 7 M (x) and satisfy the anti- 
commutation rule 7^ + 7i,7 M = 2g M ,,. The curvature independent ordinary Dirac 
matrices are denoted by 7^, where /j,,a = 0, 1,2,3; 7a7j> + 7&7 fl = 2^; T]QQ = I, 
rju = 1, rjQi = 0. The relation between 7 M (x) and 7 is given by vierbein-components 
L defined through the relation 



For R-W spacetime the vierbein are fixed through 

Here z = t, x l ~ r, x 1 = 0, jr 1 = ip and ds 2 is given by (2). We find 
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(14) 



In case of flat spacetime k (i.e. p = 1), L a should reduce to a diagonal form. But 
equation (4) does not reduce to diagonal form for k = as is expected. This is due to the 
fact that the vierbein is defined up to a Lorentz transformation in the flat tangent space. 
As Lorentz transformation for the 3-space dimensions can be viewed as a rotation, we 
rotate the vierbein (14) by multiplying it with the rotation matrix 



r 



/I \ 

sin 6 cosy sin 9 sine/? cos# 

cos 6 cos (p cos 9 simp sin# 

sin (p cos y J 



(15) 



or 



to obtain the new vierbein as, 




(15a) 



(16) 



For k = i.e. p = 1, this matrix turns out to be a diagonal one and moreover, it is also 
symmetric. 

Thus we can calculate the covariant curvature dependent Dirac matrices 7 ;i (x) in terms 
of curvature independent ordinary Dirac matrices 7^ by using equations (11) and (16) 
which are as follows: 



To (*) =70, 
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relations ^(x) = g^jvfa) and equations (6) and (17), are obtained as 



2 



- Phi - 



2.4 5pm connections 

The spin connections are given by the equation 

r / i = -i[7 I/ W,7,W||J, 

where the covariant derivative is given by 

7^0)n M = 7^W| M - r^7 Q 0). 

We find T M by using equations (7), (8), (9), (10), (17) and (18) as 
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From equations (18) and (21), we get 
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[so using equations (3) and (18) we may write 

/ x\ %2 xj\ 

p 71 + 72 + 73 }o r 

V r r r/ 
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finally, Dirac equation (1) takes the form through (22) and (23) as 
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/e a more compact form of equation (24), we use the transformation of the curvature 
endent Dirac matrices in spherical polar coordinates by 



cos#\ /7i\ 

76) = I cos cos <> cos 9 sine/? sin# 72 h 
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(25) 



^-A-(l~p)^-m^ = Q. 

(26) 

ultiplying equation (26) by 170 and using the relation 7,- = 700;; , a/ = 707;, it 
nes 



(27) 



(28) 



i reduces to the flat spacetime equation for p = 1 i.e., 
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3. Solution of the Dirac equation 



It is better to find out a diagonal tetrad gauge than the cartesian tetrad equation, becai 
is not so simple to find out the solution of eq. (27) for p 1, as the flat space 
equation for p = 0. On the other hand there is a study on the necessary and suffii 
condition for the separability of equations obtained from the diagonal tetrad ga 
Moreover, there is a unitary transformation that maps the Dirac equation in the diag 
tetrad gauge into the equation in cartesian tetrad gauge. 

In this case, metric (2) can be written by changing the variables 

dt = exp(5a(r))dr, 
a(t) =exp(ia(r)), 



as, 



such that 



=< X, for* = 0, 

[sinh%, for k = 1. 

Using metric (30), the Dirac equation (1) takes the form 



= e Q(r) [dr 2 - {dr + 2 (x)(d6> 2 + sin 2 




or simply 






where d = da/dr and = d/dx- The unitary transformation 5" that maps equation ( 
into equation (27) is 



where 



= o (^72 + 7273 + 737i + 1) exp ( - -7371 ) exp (- ^7172) , 
^ y ji y \ 2, J 



and its inverse is 



= exp (-71 72) exp (-7371 ) -(1 - 7172 -7273 - 737i)- 
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transformation some of the identities (see Appendix A in which one of the identities is 
given in detail) are given below: 

'\ / f' 

d x +\ (Sip) = (71 sin#cos(/?+72sm$sin</?+73Cos6>)lc>x+ 

S~ 1 72(<90+ \ cot 9] (Si/j) = I (71 sin Q cos yp+72 sin 6 sin<p+73 cos 6}ip 

+ ^ cot (9(7! cos 9 cos 1/5 + 72 cos sin (p 73 sin 0)ijj (36) 

+ (71 cos 9 cos (p + 72 cos sin ip 73 sin 9}d$il), 

S"~ 1 73<9 ,(5''0) = 5 (71 cos (/P + 72 sin </?)-0 + (72 cos 1/9 71 sin i 

From equations (29) and (31), we get 



~>~r . ^- 

Using these results and the identities (36), it is possible to show that equation (32a) is 
exactly mapped into equation (27) (see Appendix B) 
In order to separate the variables of (32), we take 

^ = (X) s i n 9& ' T $, (37) 

with $ = 5(7071 + 727o + 7o7s - 7o7i727s) and 

\ 

(38) 



where n takes the values 1/2, 3/2, 5/2, ____ We get the following four equations: 
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First term of equations (39), (40), (41), (42) are dependent on r only and rest parts of 
each of the equations are independent of r. Hence the first term of each of the equations 
should be equal to a constant A (say) and thus we get, for the functions/(r) and g(r), the 
coupled equations: 



L 



r) g( T } = -i\f( T } 



(43) 
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~ + imeWw /(r) - -zAg(r). (44) 



and 



Premultiplying (43) by ((d/drHime^aM) and equation (44) by ((d/dr) - / 
and using eqs (43) and (44), results in a pair of uncoupled second order equations 



and 



+ A 2 + m 2 ew /(r) = (45) 

dr 2 2 



^M + A 2 + mV^j *(r) = 0. (46) 

J 



2 



The result obtained in this section i.e. equations (43) to (46) differs from Barut and Singh 
but is consistent with Sahni. We will consider (43) and (44) to understand particle 
production in Kasner spacetime and in some other cases. 

4. Parametric resonance particle production 

We write eqs (43) to (46) in the form 

,f = H*. (47) 

and in the second order form 

[ + (A 2 + m 2 c 2 (r) imc(T)) j <0 = 0. (48) 

Here 

and H = cr 3 rac(r) + ^A, (49) 



and CT* are the usual Pauli 2x2 matrices. The complex time WKB (CWKB) approxima- 
tion allows us to calculate pair production amplitude as follows [9-12]. We determine the 
complex turning points r\^ from 

"(7-1,2) = [A 2 + m 2 c 2 (r)imc(r)] 1/2 = 0, (50) 

and the pair production amplitude is given by 

fexp[2/f'a; A (T)dr] 
- 



In obtaining eq. (51) we assumed that particle-antiparticle states are determined through 
zeroeth order WKB approximation to eq. (48). From eq. (51) it follows that R has poles at 

# = 0,1,2,.... (52) 
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When the condition (52) is violated the poles become resonances and we get resonance 
particle production. Because of recent interest in parametric resonance particle produc- 
tion [23-26], we consider Mathieu-type equation which is generally used to understand 
catastrophic particle production 

~ + [h - 20cos(2z)]w = 0. (53) 

The advantage of our method lies in the WKB vacuum definition (usually taken at 
\r\ > oo). In cosmology we sometime require to know the particle production at an 
intermediate stage to obtain the late time evolution of the universe. In such cases the 
parameter 6 may be large and the WKB vacuum definition remains valid even at finite r. 

To demonstrate, we consider large 9 approximation and replace the periodic barriers by 
piece-wise inverted harmonic oscillator; and take z\ < f such that cos(2z) (1 2z 2 ). 
We write eq.(53) as 

A2 ~ ,1. 

. + z' 2 ] = 0, (54) 



dz' 2 
with 



The turning points are at z' i\ . We evaluate eq. (51) to get 
, D ,2 exp(-TrA) 



1 + 

exp( 7rv$e) ,_^ N 

= -=. . (56) 

1 + exp( n\/9e) 

The solution of eq.(54) can be written in terms of the standard parabolic cylindrical 
function D_, i ^g e+ , l , 2 ^(e in / 4 26 l / 4 z) which in large 9 approximation will satisfy, at finite z, 
the conditions of WKB definition of particle-antiparticle states. Using the relation 
r(z)F(l z) = 7r/(sin(7iz)), we can determine the standard Bogolubov coefficients a 
and /3 from eq.(56) with identification R = (3/ot as 



(3 = Q -\ (57) 

r((i + iV&)/2) ^ . { J 

This result exactly coincides with [26]. The solution of eq.(53) is written as 

u(z) = ^P(z], (58) 

with A real and positive, and P(z), a periodic function. To obtain A in (58) we have to 
consider a large number of repeated barrier crossovers to obtain A~ (l/7r) ln(e~' rV ^ e ^ 2 ) = 
(\/^e)/2.Thus we find that parametric resonance particle production is basically CWKB 
resonance particle production in parameter space (0, h}. 

We have considered eq. (53) to understand the effectiveness of our approach. In the 
next chapter we consider the particle production in Kasner's spacetime through CWKB to 
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understand some salient features of particle production. The above example of particle 
production in parameter space allows us to define particle-antiparticle states even when 

?7| 7^ CO. 



5. The Kasner spacetime 

We now concentrate primarily on the Kasner spacetime to study the pair production. In 
this case we find that temporal equation shows the existence of complex turning points 
which are needed to evaluate the pair production. 
The Kasner metric is given by 

ds 2 = di 2 - f^'ck 2 - t^dy 2 - f^dz 2 , (59) 

where a\ , ai and 03 are real numbers. The spacetime (59) is a model of anisotropically 
expanding universe and is a solution of vacuum Einstein equation subject to the constraint 

a\ + a 2 + a 3 = a\ + a\ + al = 1. (60) 



In the non-degenerate case the <z,-s are non-vanishing and none of the a/s are equal. On the 
other hand in the degenerate case one of the a,-s is equal to unity, while the others are 
vanishing. We choose corresponding to eq. (59) 

7(0=7o=A Y(0=fV', *=1,2,3. (61) 

Also we have 

%i , = diag(l, -1,-1,-1) (62) 

and 

o ( l \ i a t i (Q a*\ i>\ 

7 =(Q ij ; ^^ cf =( a > oj' (63) 

represent flat 7-matrices, while cr' are 2 x 2 Pauli matrices. The Christoffel symbols F% p 
and the spin connections F M are evaluated to get 

7% = -^- (64) 

Using eq. (64) in equation (1) we have the Dirac equation in Kasner spacetime to have the 
form 

d + i(r flI V 0! + T fl2 7 2 a 2 + r fl3 7 3 d 3 ) + ^- - ml ^ = 0. (65) 

^t J 

Substituting 

i> = r 1 ^, (66) 

in eq. (65) and using eq. (63) we have 

[do + r fll a 1 di + r a *a 2 d 2 + r fl3 a?d 3 + rfni$ = 0. (67) 
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If we introduce 

1 



(68) 

and also using eq. (63) then, like eqs (43) and (44), the functions f\ and f 2 satisfy the 
coupled equations 

(do + im)/i + i(r ai (r l ki + r a ^k 2 + r a3 a 3 fc 3 )/2 = 0, (69) 

(do - im] f 2 + i(t-* a l k, + r a *v 2 k 2 + r<Vfc 3 ) /i = 0. (70) 

We note that we can go over to the flat space by identifying 

We seek to find the WKB solutions of (69) and (70) both in non-degenerate and in 
degenerate case. 

5.1 Non-degenerate case 
We use the transformations 

where 

_ -a, 1 - 2 2 - fla 3 ,_ /A 



A = y'c? + c\ + c\, a = r a 'ki, U = 
and 



As a result 9,' takes the form 

9 _ d ic l c 2 c 3 (a l -a 2 ) . 
d ' ^ 2A(cf + ci) r ' 

Equations (69) and (70) becomes 

($/+im)/+M / s = (74) 

and 

(df-im)g + iA'f = Q. (75) 

Premultiplying eq. (74) by iA'/X 2 and operating by (cfy zm) from left and also using eq. 
(75) we get uncoupled equation for /, 

[d? 2 - d f (In A)<9,< - imdf (In A) -I- A 2 + m 2 ] / = 0. (76) 
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Similarly tor g using j(t, m), 

g(t,m)=f(t t -m). (77) 

Substituting/ = e (1/2) / 9 <' (lnA)d// ; y in eq. (76) and also using eq. (73) it turns to be of the 
form 

[fl? + n 2 (01z = 0, (78) 

where 

v __ Ji/2) f2/ 9 (0df- (} _ c\c 2 c 2 ( ai -a 2 ) i 
y - e Z) * ( - 2A(c? + ci) ' ' 

5.2 Degenerate case 

In curved spacetime the role of vacuum is not well resolved. To understand the role of 
vacuum, we concentrate on degenerate Kasner spacetime in the light of CWKB technique 
and evaluate the pair production amplitude. This is done to elucidate the vacuum with 
respect to strong and weak vacuum [8, 27-29], The study of particle production in the 
degenerate Kasner spacetime is interesting because of its striking similarity with the 
Milne universe. 

In the degenerate case we take a\ = 03 = and a 2 = 1, so that Kasner metric (59) 
reduces to 

ds 2 = df- (dx 2 + t 2 dy 2 + dz 2 }. (79) 

Accordingly eq. (67) becomes 

[do + a l di + r l a 2 d2 + a 3 ^ + ry m]^ = 0. (80) 

We insert in eq. (80) the following 

^~e ik ' x x, (81) 

where 

X=g) (82) 

and define 



(83) 
+ (3m. 
Equation (80) then takes the form 

[ido - a 2 k 2 r l - j3'M\x = 0. (84) 

As eq. (84) has only two anticommuting variables we may choose 

/5' = cr 3 and a 2 ^ 1 , (85) 
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as if (84) reduces to Dirac equation 

idQijj = H-0, (86) 

where 

/,,. \ 

(87) 

TT IT, ^-1 _j_ f3\/i fR8l 

This equation has a remarkable similarity with Dirac equation in time dependent 
electromagnetic field [18]. The basic feature of reflection and particle production in the 
work of Cornwall and Tiktopoulos [18] is that the external field (in our case it is the 
gravitational field) acts between t\ and t 2 . For t < t\ there is no particle but at t > t 2 (or 
t > +00) the solutions behave such that it becomes a mixture of particle-antiparticle 
solutions and is written in terms of free orthonormals. Cornwall and Tiktopoulos 
demanded that between t \ and t 2 the spinor suffers a rotation of odd multiple of ?r in order 
to have pair production. It may be mentioned that the particle production in an electric 
and gravitational field was also explicitly calculated for scalar and spinors by Sahni [30]; 
however, our approach is very much similar to particle production in flat spacetime using 
Feynman-Stuckleburg prescription. 

We generalize the above approach through our complex time multiple reflection 
technique for which it is necessary to find out two solutions that will behave like particle 
and antiparticle at t oo in order to fix up the vacuum. Though not mentioned in the 
work of Cornwall and Tiktopoulos we observe that t t\ and t t 2 must be two turning 
points to be determined from the second order Dirac equation corresponding to eq. (84). 
We show that for a system having turning points this particle-antiparticle transformation 
takes place very naturally [15]. 

(i) To generalize the work of Cornwall and Tiktopoulos [18] we used particle- 
antiparticle definitions in curved spacetime as follows. We adopt WKB definition of 
particle-antiparticle [8, 27-29] solutions such that 

lim C/|j n ut ~ exp(z,$o n ut ) (89) 

and 



respectively. Here 

51 = ^5 (91) 

and S is a solution of Hamilton-Jacobi equation. 

(ii) To have reflection in time we must have turning points. 

(iii) The pair production amplitude depends on two fixed points i.e. the turning points. 

We use eq. (82) in eq. (84) to have coupled equations 

(do + iM] /i + ik 2 r l f 2 = 0, (92) 

f l =Q. (93) 
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Premultiplying eq. (92) by / and then operating by (do - iM) from left ana using (yj) we 
have, like eq. (78), 

= 0, (94) 



where we have introduced 

f = e* (95) 

and 

0(77) = [M 2 e 2? > + iMe* 1 + f.] 1/2 . (96) 

We now find the particle production amplitude |/?| 2 by our non-perturbative method of 
WKB approximation in complex time. Equation (94) is of the same form of a 
Schrodinger-like equation (not in space) in conformal time coordinate viz., 

= 0, (97) 



dif 
where 

fifo) = (or - V) 2 , (98) 

with u playing the role of energy and V being a potential term. We define the following 

r = Mt\ 



c = ik 2 . 
The turning points r\ and TI are now obtained from 

a(n, 2 ) = o, (100) 

to be 

TI = b + \/b 2 + c 1 . 

_ (101) 

T2 = b ~ \/b 2 + C 1 . 

Since the particle production amplitude R (see Appendix C) is given by 



where from equation 

%)= / > 0(7 ? )dr ?) (103) 

/ 

we get 

_, . 7T 7T^ 2 

5(ri)=-- (104) 
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S(Ti,T 2 }=-7rik 2 . (105) 

The pair production amplitude can now be obtained employing (102) to be 



2 



The result shows that there is particle production. If we wish to obtain the result for Milne 
universe herefrom, we may set k\ = kj = which means from eq. (83) 

M = m, (107) 

and the pair production amplitude will be of the same form of eq. (106), as it is 
independent of M = m and thus we get particle production even in Milne universe in two 
dimension. It is a new result and needs careful analysis. We have based our calculation on 
WKB approximation and adopted the WKB prescription of particle-antiparticle. So the in 
vacuum and out vacuum are adiabatic vacuum. It is surprising to find pair production 
even in adiabatic vacuum which is possible when |0, in) ^ |0, out) though both are 
defined in terms of adiabatic modes. The standard result for the Milne universe in two 
dimension [9, 16, 31] are as follows: 

(1) There is no pair production in adiabatic vacuum. 

(2) There is pair production in conformal vacuum. 

These two facts are based upon the solutions of the temporal equation [16]. However, 
our result does not rely upon a particular solution and depends only on the turning points 
and the extra phase obtained through reflections. It is a global feature of Robertson- 
Walker spacetime. 

The standard arguments for the Milne or Chitre Hartle universe both in two and four 
dimensions reveal that under coordinate transformation the metric reduces to Minkowski 
spacetime and hence we do not see any pair production in Milne universe. From our 
approach it is clear that as soon as one does the coordinate transformation to reduce to a 
Minkowski spacetime-like metric, (a) one squeezes the spacetime instead of the whole as 
in the original universe, (b) It also gets rid of the turning points. This in turn means that an 
electron in a Milne universe is definitely not equivalent to some other electron in another 
spacetime obtained after coordinate transformation. Hence seeing the similarity of Milne 
universe with Minkowski spacetime it is not desirable to comment on the existence or 
non-existence of pair production. 

6. Discussion 

The method which we have developed has its basis in the Feynman-Stuckleburg 
prescription according to which negative energy particle solutions propagating backward 
in time is considered equivalent to positive energy antiparticle solutions propagating 
forward in time. Accordingly the e + e~ pair creation from a potential (where we have 
positive energy e + and e~, both propagating forward in time) may be treated as being 
equivalent to the reflection of a negative energy electron initially propagating backward 
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in time, to suffer reflection in time thereby giving way to a positive energy electron 
moving forward in time. Alternatively, particle production can be viewed upon as a 
process of reflection in time due to which there is no particle present at t > oo but at 
t +00 there is a particle moving forward in time and an antiparticle moving backward 
in time. In CWKB the wave starts from a real point and suffers reflection at the complex 
turning points. Knoll and Schaeffer [32,33], Schrempp and Schrempp [34] developed the 
method of complex semiclassical paths in space which they have utilized in calculating 
scattering amplitude (in case of reflection in space). We extended their method invoking 
complex semiclassical paths in time which is used to calculate pair production amplitude 
(in the case of reflection in time). 

Elsewhere [10] we have established an analogy of CWKB with moving mirror example 
in curved spacetime. Similarity of pair production in gravitational field with that of 
electromagnetic field has been elucidated in our earlier works [12]. As the production of 
particles carried out in our approach is remarkably similar to pair production in 
electromagnetic field, it is thus supposed that instability of the vacuum is the root cause of 
pair production in such a situation as in electromagnetic case. Im ef f , which is needed to 
measure the instability of vacuum, can be calculated, as in our earlier work [10], in 
CWKB, without taking recourse to Bogolubov coefficient. 

Let us now present a synopsis of the work done here, (i) The particle-antiparticle 
model in curved spacetime is fixed by zeroeth order WKB solution of the field equation. 
In other words we have invoked Parker's definition, (ii) The vacuum is chosen from the 
CWKB solutions which serve as mode solutions for proper adiabatic vacuum, (iii) 
Occurrence of turning points ensures particle production. This fact is hitherto unnoticed 
in the literature. These turning points are the characteristic of a given spacetime. (iv) 
Occurrence of turning points is a common feature of all standard expanding spacetimes. 
Hence particle production is a general characteristic of expanding spacetimes. (v) The 
turning points are determined utilizing eq. (98). (vi) The in vacuum in our degenerate 
Kasner case is Hartle Hawking vacuum as in collapsing star problem [10,31]. The 
out vacuum is such that as the particle evolves from the in vacuum, it does not find 
the out vacuum as a physical vacuum i.e. |0, in) ^ |0, out). This occurs because 
jQ(r)dr| _ ^ 0. (vii) Mixing of positive and negative frequency modes result in 
particle production. The situation is like the blackhole problem [10, 30] where mixing of 
positive and negative frequency states occur at the horizon, (viii) Vacuum acts like a 
blackhole. In blackhole case hole mass decreases but here vacuum energy changes. In 
both cases Im c ff remains finite, (ix) The standard explanation of no particle production 
in Milne universe is that under coordinate transformation, (for two dimension) 

v = e^coshjc, 

(108) 
v = e^sinhjc, 

with < y < oo and oo < v 1 < oo, it reduces to Minkowski spacetime. Hence we do 
not expect particle production in two dimensional Milne universe. Our calculation, 
however, contradicts the situation both in two as well as in four dimension. The coordi- 
nate transformation washes out the turning points and hence there is no particle produc- 
tion. But the turning points are crucial and acts like horizon [10, 31, 35] and is responsible 
for mixing of positive and negative frequency states resulting in pair production. 



Details in this respect are discussed in our earlier works [10, 1 1]. The basic calculation is 
that the vacuum in Milne universe is not the same as Minkowski vacuum and in a 
gravitational field, the instability of this vacuum results in particle production. In complex 
time WKB approximation (CWKB) technique, the aforementioned ideas find an elegant 
manifestation. 

We have also studied catastrophic particle production under periodic perturbation 
through CWKB and find results consistent with other works. This allows us to study 
particle production in the light of CWKB to understand some aspects of cosmological 
evolution of the universe. 



Appendix A 

Among the three only one of the identities 

-r l Ts(ft + 5cot0)(styO 

= 2 (7 1 S ^ n ^ cos 'P + 72 sin 9 sin (p + 73 cos d)ijj 
-\-\ cot #(71 cos 9 cos (p + 72 cos 6 sin (^73 sin#)-0 
+ (71 cos 6 cos < + 72 cos 9 sin (p 73 sin 9) cV0, ( Al ) 

is shown in detail here and others follow the same rules. 
Here we strictly follow 



and 



S = i( 7l72 + 7273 + 737! 



Then the left hand side of equation (Al) may be written as 

S- l n(do+{cotO)(S^ = (S- l i2Sd 8 + S- } j2(deS)+cot0S- l j2S}iJ. (A2) 
Here 

STV = le^W^ (1-7 172 - 7273-7371)72 
x (7172 + 7273+7371 + 

= 1^/2)7,7^/2)737, (J _ 
X 



-7i7273-72+7i -73+73+71 +7 2 




e (>/2)7,72 
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cos 9 - 73 sin 6} (cos | - 7, 72 sin ^) 

^ ^- ^-r / 

= (cos - + 71 72 sin -J | (71 cos - 73 sin 0)cos - 

+ (72 cos 4- 71 7273 sin 0)sin -z \ 
= 71 cos 9 cos c^ 73 sin 9 + 72 cos Q sin <>, (A3) 



S~ ' 72 (505) = - S~ 1 72 (71 72 +7273 +737i + 1 )00 cos - - 7371 sin - 

= 16^/2)7.7^(0/2)737, (1 _ 

1/91 

~ 2 Sin 2 ~ 2 7371 2 



cos 
8 2 



+7371 sin? VTI sin?+7 3 cos? )e~ (v3/2)7172 
2J Y 2 2y 

1 72 (73 cos ^ + 71 sin 9) (cos ^ - 71 72 sin -) 

\ jL / 



_ _ 

j-f 

~ ~ 2 ( C S f + 7l 72 Sin f ) i/ 73 COS ^ + 73 Sin ^ COS f 
+(-717273 cos 9 + 72 sin 9) sin - j 

= ~"2 (73 cos + 7i sin ^ cos <yO + 72 sin sin </?) , (A4) 

1 cot ^5~' 72^=5 cot ^(71 cos<9cos<p 73 sin^ + 7 2 cos^sin(/?). (A5) 

Using equations (A3), (A4) and (A5) in equation (A2) give the required result of 
equation (36). 

Appendix B 

The unitary transformation (34) maps equation (32) into (28) which has to be proved 
here. 

From equations (29) and (31) we get 



v / sj \ / i i A * 

Using these with a,- = 707, and the identities of equation (36), equation (32a) becomes 

/ 3 -\ f 

e (i/2) Q (r) (Q ( + -- J -7o< (71 sin B cos (+72 sin 9 sin ^+73 cos 5) (pd r + -) 

\ jw Ci- / I \ / / 

(71 sin cos <> + 72 sin 9 sin ^? + 73 cos 9} 

1 / cos 2 B cos 2 . 






rsin0 



(-71 sin <y9 + 72 cos (p)d v > + zm7 e (1//2)Q(r) 



or 



d t -\ < (71 sin 9 cos <p + 72 sin 9 sin (p + 73 cos 6) (pd r + 

2 a a I \ 



(71 sin cos </? -f 72 sin sin (p + 73 cos 0) 



+ - (71 cos cos (p + 72 cos 6 simp 73 sin 9}d0 
r 



1 



r sm 



(~7i sin ^ + 72 cos 



This equation becomes, after using the transformation (26), 

a, + + to - 7, (^ + ^) - ;7r + \1ide 



and thus the result. 

Appendix C 

We determine the turning points, using eq. (98), from 



- 0. 



Defining 



and 



the boundary condition takes the form 



f] > -CO 
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where R is the reflection amplitude. In eqs (C4) and (C5) 77 is real but the turning points 
are complex. In CWKB the reflection amplitude is identified as the pair production 
amplitude. The CWKB solution of eq. (97) is 

V^ r - e iS(T>i ,Ji)l 2 M (C6) 



where 77 is real and 770 is real, arbitrary and 770 > f\- In (C6), the first term represent a wave 
starting at 770 and moving left reaches real 97 < 770. The second term takes the contribution 
of the reflected part as follows. A wave starting at 770 reaches the complex turning point 771 
and after bouncing back from 771 reaches 77. This is represented as 

(-Oexp[zS(T7i,7)d)-tfMi)]. (C7) 

This contribution (C7) to the second term of (C6) is then multiplied by the repeated 
reflections between 771 and 772 and so comes the following factor in the second term of 
(C6), namely 



1 IT- (C8) 



n=o 



For convenience we have neglected the WKB pre-exponential all through. Comparing 
(C6) with the boundary condition (C5) we have 



If we consider exp[z,S(77, 770)] as antiparticle solution then (C6) is interpreted as the 
reflection of antiparticle from the turning point 771 and the reflected part is interpreted as 
particle moving forward in time. This is Klein paradox-like situation not in space but in 
time. 
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1. Introduction 

Study of symmetry changing phase transitions of a quantum field in the presence of a 
surrounding thermal bath is very important in the study of the evolution of the universe 
and in the analysis of very high energy collisions where very high matter and radiation 
density exist. Detailed study of these phase transitions has been done by various 
authors [1-4]. The effective potential method is very useful in studying spontaneous 
symmetry breaking (SSB) at zero temperature [5,6]. Estimation of the critical 
temperature of phase transitions can be done by extending this approach to finite 
temperature. These studies mainly using loop expansion techniques have played a pivotal 
role in framing our understanding about the early universe, unified theories, quark gluon 
plasma etc. 

Recently there has been a revival of interest in finite temperature quantum field theory, 
apparently caused by the recognition of the importance of symmetry breaking phase 
transitions and the problem of precise determination of critical temperature. An accurate 
analysis of phase transitions (both analytical and numerical) becomes necessary because 
most of the cosmological models critically depend on it. For example baryon asymmetry 
may be generated at the electro- weak level if the phase transition is of first order [7-10]. 
For calculation of critical temperature various perturbative and non-perturbative 
techniques have been suggested [11-19]. 

Effective potential defined as single Legendre transform provides an efficient way 
to obtain quantum corrections to the classical potential. But this popular method 
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suffers from a serious shortcoming. Since it is defined as a single Legendre trans- 
form, it is always a convex function. This forbids the double well shape for the 
exact effective potential and implies the absence of local maximum at the symmetric 
origin. But in a theory possessing SSB classical analysis predicts a maximum at 
the origin. Various procedures to avoid this difficulty have been suggested earlier 
[3,20,21]. One of the most efficient of these approaches is to define an effective 
action by including a source K(x,y) coupled to a term which is quadratic in the 
field variable. By this procedure an effective potential with a proper loop expansion 
to each order which is not convex, is obtained. This idea was first put forward 
by Hawking and Moss in the context of quantum field theory in the early universe 
[22]. A self-consistent improvement for the finite temperature effective potential has 
also been suggested [23]. In this formalism it is possible to sum a large class of 
ordinary perturbation theory diagrams that contribute to the effective action and the 
gap equation which determines the form of the propagator, is obtained by variational 
method. 

All the above improvements developed for conventional finite temperature effective 
potential are based on an important contribution by Cornwall, Jackiw and Tomboulis 
(CJT). They defined the effective action for composite fields in flat space and zero 
temperature as a double Legendre transform with two sources J(x] and K(x, y). These 
two sources are coupled respectively to 0(jc) and (j>(y) [24]. The CJT formalism is 
considered to be best suited for studying phase transitions because it uses a generalized 
effective action in which not only the mean field but also the correlation functions appear 
as independent variables. A simple series expansion has been developed for this improved 
effective action [24]. 

The CJT formalism has recently been used to resolve various difficulties in quantum 
field theory [25, 26]. For example, it has been applied to the triviality problem in $ 4 
theory [27]. An improved effective potential based on this formalism has also been 
developed [28]. The self-consistent improvement of finite temperature effective potential 
involves the summation of daisy and superdaisy diagrams, for which a novel re- 
summation procedure has been recently proposed [29]. 

Recently, there has been considerable interest in the application of functional 
techniques in field theories with dimensions less than (3 + 1), mainly because some of 
the problems afflicting 4-dimensional theories are absent there [30-32]. According to 
Coleman's theorem, spontaneous symmetry breaking can occur only when the 
dimensions are higher than (1 + 1). In (2+1) dimensions the most general 
renormalizable theory is for a $ 6 model. This model has been studied earlier by 
using various methods and it has been shown that it possesses an ultraviolet fixed 
point in I/TV expansion and Gaussian approximation [33-35]. Finite temperature 
field theory has also been analysed [17, 18]. It is well known that in finite temperature 
CJT analysis of <fr 4 theory, the effective potential shows a cut-off dependence due to 
the presence of a (X(j) 4 /l2) term. It is natural to think that in lower dimensions 
where coupling constant renorrnalization is not required this difficulty will not be 
there. Even though this is found to be true for <fr 4 theory in (2+1) dimensions, an 
unrenonnalized mass term appears in the expression for effective potential for <fr 6 
theory. Thus the difficulty persists in a disguised form, in (2+1) dimensional $ 6 
theory. 
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2. CJT formalism 

The CJT method provides a generalization of the ordinary effective action F(0) (the 
generating functional for single particle irreducible n-point functions). This generalized 
effective action F(</>, G) depends both on <j)(x} the expectation value of quantum field 
$(*), and G(x,y] the expectation value of time ordered product T(3>(x)$(y}}. Physical 
(on-shell) solutions require the following variational equations: 



(2) 



6G(x,y) 
Consider the vacuum persistence amplitude Z(J, K] in the presence of two source terms 



Z(J,K] =N 

J *" 

(3) 
W(J,K), the generating functional for connected diagrams is defined as 

Z(J,K)~exp[iW(J,K)}. (4) 

The classical action /($) = / d*x(x) may be written as 

(X)DQ ' (x - y)0(y) + / int (0) (5) 



(6) 

y) is the free propagator that satisfies 

/>o 1 (*-y) = -(D+2) 4 (*-y). (7) 

The generalized effective action r(0, G) is the double Legendre transform of W(J, K}, 

, G) =W(J,K) - I d 4 W*(*)y(jc) - \l d 4 xd* y [<S>(x}K(x,y)3>(y}} 

y), (8) 



(9) 



where J(x) and K(x,y] are determined by 

SW(J,K) _ ^ 



6J(x} 

-G(x,y)]. (10) 
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By actually performing functional differentiation on (8) we find 

= _y W _/>,*(,,,)#>,), (ii) 

x) J 

, G) l ,\ 



In the absence of sources, (1) and (2) are regained which permit a variational solution. 
The conventional effective action F(</>) = F(0, GO) where GO is the solution of (2). 
Generalized effective action F(0, G) is the generating functional for the two particle 
irreducible (2PI) Greens functions expressed in terms of the full propagator. The series 
expansion for F(0, G) is shown to be [24] 

) =7 c]ass .(0)+ itrlnDoG- 1 + ^tr[D~ l G - 1] + F^(0,G), (13) 



where 'tr' is the functional trace 'In' is the functional logarithm and D~ { G is the 
functional product. The inverse propagator is defined as the functional second derivative 
of action: 



Computation of the quantity F^ 2 ^ (0, G) is done as follows. In the classical action /($), the 
field $ is shifted by 0(x). The shifted action /($ + </>), possesses terms cubic and higher 
in $. These define as interaction part /; nt with vertices depending on 4>(x). T^(cj)^ G) is 
then given by all the 2PI vacuum graphs and the propagator is set equal to G(JC,J). The 
theory in its full generality is not translationally invariant since vertices depend on <f>(x) 
and G is not a function of (x y) alone. The propagator of the theory is determined by 
finding the gap equation for G using the variation equations [eqs (1) and (2)]. 

3. Effective potential 

Translationally invariant solutions are obtained by imposing the following conditions 
(homogeneous states): 

(j>(x) - Constant, (15) 

G(x,y}=G( X -y), (16) 

(17) 

where (0, G) is the minimum of the energy when varying over all the normalized states 
with constraints: 

<*(*)) = #*), (18) 

(x,y) t (19) 

- l x, (20) 



where v is the space-time dimension of the theory. Thus the effective potential is 
given by 

, G) 

(21) 



^ . 

A series expansion for the effective potential is obtained by defining the following 
Fourier transformed propagators 

= / d*jce*<*-tfG(jc - y), (22) 

(23) 

(* - y), (24) 

, G) = l/(0) +^/~ In detDoWG- 1 ^) 

/:-l] + V (2) (^G), (25) 



where /(</>) is the classical potential, V( 2 ) (</>, G) is the sum of 2PI vacuum graphs, with 
vertices given by /j nt (0, $) and the propagator is set equal to G(k). The field on which 
vertices depend is now a constant parameter. Trace and logarithms apply to component 
fields and determinants are no more functional. 

To describe the theory at finite temperature we use the Euclidean time r satisfying the 
boundary conditions < r < (3 = l/T. All the Feynman diagrams (2PI diagrams) 
developed at zero temperature are valid here also. The Feynman rules for writing the 
algebra of the diagrams are different at finite temperature [21]. They are 

(26) 
loop integral - --, (27) 




vertex delta function - (3 (2ir) "~ ' 8 ^ <5"~ l [ ) J k t ] . (28) 

u n 

Field $ satisfies the periodic boundary conditions 

(-!-) "(14 

With these modifications we can write the series expansion for finite temperature CJT 
effective potential (analog of Gibb's potential) with time integration suppressed and a 
summation performed 

1 r -iv-lr, 

h^r- (3) 



Pramana - }. Phys., Vol. 50, No. 2, February 1998 137 



K P Satheesh and K Babu Joseph 



4. 0> 6 Theory 

The classical potential of the theory is given by 

r r / -n\ 1 1^1 A]J - 4 CB -r fi / o 1 \ 

C/Y <P ) = m^. Q? -) y? + - - $ ( J I ) 

where A and represent coupling constants for <$> 4 and $ 6 interactions respectively and 
the subscript B indicates bare parameters. The functional operator D~ l is given by 

J (x-y}. (32) 

After shifting the field 3? to (3> + 0), the interaction Lagrangian takes the form 

f C "I 



A few of the 2PI vacuum graphs up to three loops are shown in figure 1 . We select only 
those graphs with vertices depending on first order in the coupling constant. This 
approximation (Hartree-Fock) corresponds to a systematic variational procedure and is 
superior to commonly used one-loop approximation. No graphs with internal lines appear. 
Daisy and superdaisy graphs shown in figure 3 are of this type. Thus only the graphs 
shown in figure 2 need be summed. The sum of the relevant 2PI graphs takes the form 



xG(x t x)G(x 1 x)+ d v xG(x,x)G(x,x)G(x,x). 

6! J 

(34) 

The expression for the finite temperature nonlocal composite operator effective action in 
Hartree-Fock approximation becomes 



10 



(35) 




etc. 3 loop 



Figure 1. 2PI graphs for $ 4 and <& 6 . 
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Figure 3. Daisy and superdaisy graphs for <& 4 and <& 6 . 
By performing variation with respect to G the modified gap equation is obtained as 

( X -y}. (36) 



By iteration it can be seen that G generates all daisy and superdaisy graphs of figure 3, 
which reveals the fact that we have achieved a definite improvement in the resummation 
of diagrams. 

Since we are interested only in translation invariant theories we fix an ansatz for G and 
define the Fourier transformed propagators [eqs (22), (23) and (24)] 

(37) 

-.. (38) 



Here the propagator is chosen in terms of an effective mass M which acts as a variational 
parameter. Effective potential in terms of M 2 and </> can be written using static 
configuration and constant background field 



Putting B = we get $ 4 theory in (2+1) dimensions. Comparing the expression obtained 
with the Gaussian effective potential studies [16] we see that both are formally similar. 
But as far as $ 6 theory is concerned dissimilarities emerge because of the term <^. It is 
interesting to note that this factor is not contributed by the daisy or superdaisy diagrams, 
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but by a graph with vertex not proportional to f . This graph is not considered in Hartree- 
Fock approximation. 

Since the effective potential is an ordinary function (not a functional) stationary 
requirements with respect to 4> and M 2 is obtained by ordinary differentiation 

8V , I T AR ,-> . CB ,A . A B ^./ \ , <B x. r ). / ^,l_. _A I n fAQ\ 



= 0. (41) 



Conventional effective potential is defined at the solution of (41). The effective mass is 
given by 

r \ /- \ t- ~] 

(42) 

Required expression for the effective potential is obtained by replacing the effective mass 
M by M((j>} in (39). Equation (39) shows certain very important peculiarities of $ 6 and $ 4 
theories relevant at zero temperature 

\ i 

(43) 

For 3> 4 theory M 2 (<) is intrinsically positive. Hence if AB < only solution to the above 
equation is cj> = (or potential is unbounded from below). That is for negative AB non- 
zero turning points do not exist. In the case of <fr 6 theory 

i-~G 2 , (44) 

non-zero turning points are possible also for A 3 < 0. <fr 6 theory in Hartee-Fock 
approximation requires up to three loops for obtaining the effects of coupling. We have 
four parts for the effective potential 

= V + V 1 +V 2 + V 3 , (45) 

(46) 



(47) 

(48) 
(49) 



5. Renormalization 

The effective mass term M($) defining the effective potential is divergent mainly due to 
the presence of G(x,x). Following renormalization prescription is employed in (2+1) 
dimensions to regularize Af (0) [6, 23]. Define 



,<. n v 

vtfW^) (50) 

with 



. (51) 

In (2+1) dimensions coupling constant renormalization is not required. Define 



r 2 

r f d 2 k 1 /AN 

/i = / -r = lim - , 

1 7 (27r) 2 2fc A-KX,^,/' 

' = G(M(</>)). Using the summation procedure developed by Dolan and Jackiw [2, 23] 
the summation in time co-ordinate can be performed 

*L ! (54) 

O _-. / iQC 1 < \ " * ' 



By actual evaluation, introducing a cut-off parameter A 

/. ( 55 ) 



Equation (54) shows that G(M0)) is the finite part of the vacuum propagator. A finite 
expression for Af ($) is obtained by expressing it in terms of the renormalized parameters 



(56) 

Second derivative of the tree level potential is defined as m^ (tree level mass) 

(57) 

(58) 



At zero temperature the second term vanishes and the last term is cut-off dependent. 
Cancellation of this divergence is obtained by combining V, V 2 and V 3 : 



(59) 
2 
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with 



(60) 



= reproduces the result of Camellia and Pi. Using unrenormalized gap equation we 
combine V, V 2 and V 2 and writing them in terms of renormalized parameters 



(61) 
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Figure 5. Effective mass for various couplings-^ 4 theory. 



In the case of (2 + 1) dimensional <]? 4 theory F(0) = A/12 which is finite. Thus unlike 
(3 + 1) dimensional <& 4 theory the effective potential does not contain any 
unrenormalized parameters. But in the case of $ 6 theory F(0) contains m which is an 
unrenormalized mass parameter. But here we can make F(4>] = by adjusting the 
parameters suitably and make the unrenormalized parameters vanish. 



6. High temperature expansion 

Evaluation of the effective potential at high temperature up to one loop level has been 
done earlier in (3 + 1) and (1 + 1) dimensions [2, 35]. Additional terms appearing in the 
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Figure 6. Effective potential for various couplings-^ 6 theory. 

xpression for the effective potential can be obtained by evaluating G(M(0)), 
< 1. The relevant integral (eq. (50)) is of the form 



h n (y] 



1 



1 



m 



hese integrals satisfy the differential equation 



dy n 

ligh temperature expansion for the integral is obtained by using the identity 

_J_ == 1 _I 4 .9V^ Z 

e 2 - 1 z 2 



\2 ' 



Multiplying the integrand by a factor x 6 for convergence, performing term by 
ntegration and letting e > at the end we get 



IO 

for 

(62) 
(63) 
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Figure 7. Effective potential for various couplings-^ 4 theory. 

7 = 0.5772 . . . which is the Euler constant. 



Other values of h can be found by using the differential equation (64). We have 
1 r d 2 k 1 1 r kdk 1 



-( 

PJ 



By using the formula (67) we get 



(67) 



(68) 



(69) 

M-71 71 

which can be used to evaluate V e ff at high temperature using (60). High temperature 
expansion for the effective mass can be obtained from (56). 

7. Conclusions 

A self-consistent improvement for the finite temperature $ 6 theory is obtained as an 
extension of CJT formalism. Certain peculiarities of the <& 6 field theory in (24-1) 
dimensions are analysed. <E> 4 theory in (2+1) dimensions does not contain any 
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Figure 8. Effective potential at various temperatures-^ 6 theory. 

unrenormalized terms unlike its (3 + 1) dimensional counterpart. < 6 theory in (2+1) 
dimensions contains divergent terms in the form of an unrenormalized mass parameter in 
the expression for effective potential, but can be made to vanish. In this model physically 
meaningful stable theory is possible both for positive and negative A, indicating the 
possibility of bound states. High temperature expansion for the effective potential is 
obtained. 

Behavior of the effective mass can be clearly understood by numerically solving the 
equation for effective mass. Figures 4 and 5 show these graphs for certain relevant values 
at zero temperature for <& 4 and <& 6 theories. The straight region parallel to the axis 
indicates imaginary values of the effective mass. It is clear that this region indicates a 
broken symmetry phase. Comparison between <I> 4 and <[> 6 theories shows that they are 
identical in shape except for higher numerical values for $ 6 model. A study of the 
behaviour of effective potential (figures 6 and 7) for different couplings is also given (in 
weak coupling range since the reliability of the approximation in the strong coupling 
range is not well established). Shape of the graph is same for a reasonable range of 
couplings with a notable difference in the stretch of the straight region. 

Graph for effective potential at finite temperature (figures 8 and 9) shows that as 
temperature increases the minimum gradually dissappears indicating an approach to 
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Figure 9. Effective potential at various temperatures-^ 6 theory high temperature 
expansion. 



Using the numerically obtained value of the effective mass, effective potential can be 
calculated using high temperature expansion (figure 9). The graphs clearly indicate 
broken symmetry phase at zero temperature. Comparison with the graphs obtained earlier 
[23] shows that behaviours of <!> 6 and 3> 4 theories are similar apart from considerable 
differences in numerical values. 
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1. Introduction 

We investigate here the effect of exchange of confined gluons among relativistically 
confined quarks in nucleon-nucleon (NN) scattering calculations. 

In our previous work [1-3], we investigated the ! 5o and 3 S\ (without S-D coupling) 
phase shifts by using different size parameters with and without the exchange of cr and IT 
mesons in the hamiltonian. The agreement of the results of both the models was 
interpreted as the 'Cheshire cat principle' [1] valid also for NN interaction. The term in 
the confined one gluon exchange potential (COGEP) arising from the confinement of 
gluons plays the role of a by giving the required attraction in the NN interaction. In our 
subsequent work [3] we studied the phase shifts of higher order partial waves without 
invoking a and TT mesons in our model with only cx s (strong coupling constant) as a new 
parameter. Even though the calculated phase shifts for higher partial waves did not 
compare well with the phenomenological phase shifts of Arndt et al [4] the proton-proton 
(pp) and neutron-proton (np) differential cross sections (&(&)} obtained were in good 
agreement with the experimental cr(G). The purpose of the present calculation is to 
establish our earlier claim of 'Cheshire cat principle' [1] of NN interaction. The Cheshire 
cat models are the version of the MIT bag model [5-7], where physics become 
completely insensitive to the change in the bag radius. In these models the shifting of the 
bag wall has no physical effect. Therefore, the radius of the bag can be pushed to infinity 
resulting in a purely fermionic description (i.e. QCD) or to zero, so that the description is 
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entirely bosonic (i.e. some effective meson theory). To test the validity of the 'Cheshire 
cat' syndrome we have performed the calculations with the exchange of a and TT mesons 
in the spirit of non-relativistic quark models (NRQM) [8-11] and with the quark core size 
0.5-0.6 fm. It is to be noted that this is smaller than the physical size of a nucleon and has 
to be supplemented by a pion cloud as in NRQM. 

For the confinement of quarks we are making use of the relativistic harmonic model 
(RHM) [12] which has been successful in explaining the properties of light hadrons. For 
the confinement of gluons we have made use of the current confinement model (CCM) 
[13-15] which was developed in the spirit of the RHM. The CCM has been quite 
successful in describing the glue-ball spectra [14]. The full hamiltonian used in our 
investigation has the usual kinetic energy, two body confinement potential and the 
COGEP along with a and TT meson potentials. In consistency with our earlier work [1] we 
have used a smaller oscillation size parameter so as to include a and TT mesons in our 
calculation. The parameters used in our calculation are listed in reference [1]. We have 
employed the exact analogous formulation of NRQM, using the resonating group method 
(ROM) [16, 17] and have calculated the scattering phase shifts as a function of energy 
using the Kohn-Hulthen-Kato variational principle formulated by Kamimura [18]. 

Here, we report the results of NN polarization observables (PO) by including phase 
shifts of all significant partial waves up to angular momentum (L) = 2 with the inclusion 
of a and TT mesons exchange potentials but taking into account the effect of confinement 
of gluons. Further, we compare the np PO calculated from the theoretical computed phase 
shifts with the existing experimental data so as to test directly the predictions of our 
calculation with the observable quantities. 



2. Spin polarization in NN scattering 

In this section we present the results of calculated PO in comparison with the PO obtained 
from the phase shifts of Arndt et al [19]. In our work we have employed the formulation 
and the notation of Bystricky et al [20]. In the formulation the polarization matrix X pq ik 
with four indices, which correspond to the polarization directions of the scatterer, recoil, 
beam and target. If the polarization of the particle is not analysed (or unpolarized) the 
corresponding index is set to zero. If the analysis of the polarization of particle is carried 
out, then these index are set to n, m or / corresponding to the spin projection on the basis 
vectors. The polarization matrix X pqi h is defined through the NN scattering matrix M. The 
scattering matrix M for each iso-spin state is expressed in terms of the invariant scattering 
amplitude a, b, c, d and e which are complex functions of energy and scattering angle 6. 
The five amplitudes a, b, c, d, e with the T-matrix (7-matrix = 5-matrix 1) elements can 
be determined from the phase shifts and the mixing parameters using the partial wave 
decomposition [20-22]. The details can be found in standard reference of Bystricky et al 
[20]. We have calculated PO in the energy range Ei ab OMeV to 250 MeV for np 
interaction. The standard convention is to call P a polarization, and D a depolarization. 
The expressions for the following variables are given in references [20-22]. 

In figure 1 we display spin observables /^ooo at 50 and 150 MeV for np interaction. The 
agreement between experimental data [19] and theoretically computed phase shifts is 
good at 150 MeV. But there is atypical deviation to the theoretical curve at 50 MeV. This 
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Figure 1. The calculated np polarization observable P,,ooo as a function of scattering 
angle (9) in degrees, in comparison with Arndt et al [19]. 



means at low energy the phase shifts are not well determined. This implies that 
polarization is sensitive to NN interaction at small distances. Similar conclusions have 
been arrived by Macharvariani et al [21]. Figure 2 is a plot of depolarization D against the 
scattering angle (0). The D determines the amount of polarization perpendicular to the 
scattering plane converted into the same after the second scattering. Here, the agreement 
between theory and experiment is better at all energies [19]. In the next section the 
important results are stated and summarized. 



3. Conclusion 

In this paper an attempt has been made to understand the NN PO in the frame work of the 
RHM and the CCM by employing the ROM formulation. The aim was to understand the 
role played by the confined gluons on NN interaction. We have computed the phase shifts 
of all partial waves up to L = 2 taking into account the confinement of gluons but by 
invoking a and TT mesons in the spirit of NRQM. With a smaller oscillator size parameter 
in consistent with RHM, the model is able to reproduce the phase shifts of most of the 
partial waves. The np PO are also in agreement with the corresponding experimental 
observables. However, for the majority of the PO no direct data exists. The confined 
gluons along with conventional a and TT mesons satisfactorily explain NN observables. 
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Figure 2. The calculated np polarization observable Z) n o n o as a function of scattering 
angle (9] in degrees, in comparison with Arndt et al [19]. 



The results obtained clearly establishes the validity of 'Cheshire cat principle' in NN 
interaction. 
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1. Introduction 

Recently the decay mode Afc-^A//'0 has been observed at the Femilab [1] and the 
production cross section times branching fraction for the decay Ab^AJ/ijj relative to 
that of B-^KJ/i/j has been measured to be 0.27 0.12 0.05. Based on the factoriza- 
tion approach, we have shown that the CDF data [1] for the ratio Br(A^ *AJ/-0)/ 
Br(fi * KJ/i/j] can be well accounted for in the framework of heavy quark effective 
theory (HQET) [2-5]. Much of the success of HQET has been in the treatment of 
decays from one heavy flavor to another, namely b c transitions. However, one 
approximation used recently in the literature has been to treat the strange quark as 
heavy, so that the decays we are interested here can be considered as heavy to heavy 
transitions. In fact the previous challenging studies [6-9] suggest that the heavy 
quark effective theory is still effective for s-system, if a constituent quark mass 
is employed for m s and l/m s corrections are taken into account properly. Therefore 
it is interesting to examine the applicability of the heavy flavor symmetry in the s- 
quark system even if the s quark mass is not large enough compared to the QCD 
scale AQCD- There is a lot of altercation in the literature regarding the fact, 'whether 
HQET could be applied to s quark system or not'. In a recent letter Chakraverty et al 
[10] have shown that heavy quark symmetry is not reliable for the s quark system. On 
the other hand applying HQET to B > K^ transitions, Robert and Ledroit [11] 
suggested that the heavy s limit may give an acceptable description of unpolarized 

155 



data but does not give reliable results when applied to polarization observables (T L /T}. In 
this paper we assume the validity of HQET for s-quark system and examine the ratio of 
the branching fractions for the decays A/, > A //?/; relative to B Q >KJ/i(j, with the 
factorization approximation. 

Factorization in two body non-lep tonic decays of heavy pseudoscalar mesons [12] was 
resurrected a few years ago as a means of estimating their decay rates using the existing 
calculation of the semileptonic decay form factors. It relates the complicated non-leptonic 
decay amplitudes to products of meson decay constant and hadronic matrix elements of 
current operators similar to the ones encountered in semileptonic decays. This 
approximation appears to work well phenomenologically in the B decays, where it has 
been tested very well [13]. 

Nonleptonic decays of heavy baryons in general receive contributions from W 
decay and W exchange diagrams. The calculation of W exchange amplitudes has to 
rely on phenomenological models and introduces large theoretical uncertainties [14]. 
The decay mode A/, AJ/'ifj is however free from such contributions and proceeds 
through the W decay diagram. This allows to make rather clean theoretical prediction for 
the decay rate. Assuming factorization, the matrix element of the effective Hamiltonian 
is given by the product of the baryonic decay transition and the meson decay 
constant, i.e., 

. (1) 



To evaluate the hadronic matrix elements we use the results of HQET, treating .s-quark as 
heavy. Since we are dealing with the evaluation of the ratio of the branching fractions for 
B +KJ/1JJ and A& >A.//^ decay processes, we have not included l/m s corrections in 
our calculations. As these small corrections appear both in numerator and in denominator, 
it is assumed to be cancelled without effecting much the leading order result. Therefore at 
the leading order these matrix elements can be described by the single universal form 
factor, the well known Isgur-Wise (IW) function. HQET predicts that the Isgur-Wise 
function is related to the overlap of the light degrees of freedom wavefunctions of the 
heavy particles. In general the light component is very complicated, but in the valence 
quark approximation the light component of a heavy meson is a light antiquark whereas 
the light component of the heavy baryon is a light diquark. This implies the fact that the 
mesonic and baryonic Isgur-Wise functions are quite different from each other and need 
to be evaluated separately. 

The IW functions are not calculable from perturbative QCD. On the contrary, heavy 
quark symmetry tells us only the normalization of the Isgur-Wise function at the zero 
recoil point i.e., when the initial and the final heavy quarks have the same velocity. Since 
apart from the zero recoil point HQET does not predict the shape of the IW function, it 
makes sense to calculate it in a model, which gives the same result in the appropriate limit 
as predicted by the HQET. Here we therefore use the quark model of Ali et al [8] to 
evaluate the mesonic Isgur-Wise function and the bound state soliton picture [15] for the 
baryonic function. 

The outline of the paper is as follows. In 2 we present the formalism for the 
description of nonleptonic decays of heavy mesons and baryons using the factorization 
approximation. The Isgur-Wise functions are calculated in 3. Section 4 contains result 
and discussions. 
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1 - 75)0 

)], (2) 

where GF is the Fermi coupling constant; Ci and Ci are the Wilson coefficients which 
contain the short distance QCD corrections. In the leading logarithmic approximations 
they are given to be [16] 

= 1.12 and C 2 (m fc ) = -0.26. (3) 



The effective Hamiltonian describes two classes of nonleptonic decays. The first class 
is of the type: A b *A<;D s and B->DD S , while the second class represent A. b -*AJ/ip 
and fi KJ/i/j; J/4> is the charmonium state. 

To evaluate the matrix element of the effective Hamiltonian, we employ factorization 
assumption. By Fierz rearrangement, we rewrite the effective Hamiltonian in a form 
suitable for the use of this assumption. Upon Fierz rearrangement, the factorized 
Hamiltonian for the decays is given as 

l - 7s )c), (4) 

where 2 = C 2 + C\/N C , with N c is the number of colors. One may expect N c =3, but 
the decay rate obtained from factorization is too small for N c = 3, while the limit N c ^oo 
gives a remarkable value. Since in this paper we are dealing with the ratio of the 
branching fractions, the dependence on N c is unimportant. 

As pointed out in the previous section we evaluate the matrix elements of the effective 
Hamiltonian using the factorization approximation. For B-^KJ/i{> transition, we there- 
fore write the transition amplitude as 



A(B(v l ] 

(5) 



The hadronic matrix element (K(^}\^(\ - 7 5 )^>i)> is extracted from the corres- 
ponding semileptonic decays of B mesons into the Kaon state. In the heavy s-lirmt it is 
given in the HQET [4] as 



where ( Wl ) is the mesonic Isgur-Wise function with <* = (v, t/,). The remaining 
matrix element can be written as 



- 75 = ^e, 
where ft is the charmonium decay constant and ej denotes its polarization. 
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Using eqs (5-7) we obtain the partial decay width ofB-*KJ/ip decay mode in the rest 
frame of initial B meson as 



(8) 



where p,j is the c.o.m momentum of the emitted particles and is given as 



IPil = K*4 -Ml- Af J) 2 - 4MJM 2 ] 1 / 2 . (9) 



Now similarly for the A*, AJ / A/J transition, the hadronic matrix element of the 
factorized amplitude from eq. (1) can be written as 

v 2 ,s))=77(^)MA (*/ 2 . ^b^O -75>A>2,.y), (10) 



where 77(0)2) is the baryonic Isgur-Wise function with u)2 = V2-v' 2 . The spinors present in 
the above equation MA(?/ O , s') and MAbC^i*) are the usual Dirac spinors for the A and A& 
states, with spin sum given as 

u(v,s)U(v, S )=^. (11) 



With eqs (1), (7) and (10) we obtain the partial decay width for A/, > A//^ decay as 



x [(Ml - Ml) 2 + Ml(Ml b + Ml - 2MJ)], (12) 

where |p 2 | is the c.o.m momentum of the emitted particles in the rest frame of A& baryon. 
After evaluating the partial decay widths for the two types of decays, we now obtain 
the ratio of their branching fractions as 

Br ( A, - A/M r Afe r(A, -> AJ/i/i] 



| P2 | 

' ^ " 



Pil 2 (wi) Mi b [(a; 2 - \}m b m s (M K +M B ) Z ] 

Thus one can see from eq. (13) that the ratio of the branching fractions is independent of 
any model dependent parameters. The only unknown quantities present in the expression 
are the Isgur-Wise functions (u>i) and 77(0^) These universal functions are evaluated in 
a simplified manner as presented in the following section. 

3. Evaluation of the Isgur-Wise functions 

The Isgur-Wise function which represents the nonperturbative QCD effects, are a 
measure of the light cloud (spectator quarks) rearrangement around the heavy quark 
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subsection 3.2. 
3.1 



The Isgur-Wise function calculated in the quark model, predicts the decay widths for rare 
B-*K*^ [8] and the branching ratios of nonleptonic B mesons [7, 17] very well. In this 
model the function (o>i) can be obtained from the overlap integral as 

(14) 

where <p\ and </>p refer to the radial wave functions of the initial and final mesons 
respectively, y'o is the spherical Bessel function of the zeroth order. The 'inertia parameter 
A' is taken to be [18] 



8 

for the constituent quark mass values m<* = 330 Me V and m s = 550 MeV. To calculate the 
Isgur-Wise function (14) we use the wave functions of harmonic oscillator for <$>\ and 0p 
in the form 

3 1/2 

exp(-/3V/2), (16) 



with oscillator strength /?. The values of /3 K and j3 B fitted in ref. [19] (J3 K = 0.34 GeV and 
{3 B = 0.41 GeV) are not equal. However we have assumed (3 to be same, for the initial and 
final mesons and obtain 



We have taken j3 = 0.295 GeV, which is extracted from the best fit in ref. [7]. 

The value of ui\ is determined by considering the kinematics of the system. Since 
we are dealing with the two body decays B(v\) K(v\}J/ip(p^), momentum conserva- 
tion gives 



(18) 
which gives 



. 
2M B M K 

Thus with eqs (17) and (19) we have obtained the mesonic Isgur-Wise function (k>i) 
to be 

0.317. (20) 
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3.2 

Here we have presented the evaluation of the Isgur-Wise function in the same manner as 
suggested in ref. [15]. The Isgur-Wise function calculated in this manner explains very 
well the two body non-leptonic decays Ab->A c P(V] [20] and the radiative rare decay 
A/, > A7 [21]. In this model the heavy baryons are treated as the bound state of the 
soliton with the heavy meson. The Isgur-Wise function is given as 

77(012) = j d 3 q$*(q)$(q + m fl (v 2 - v' 2 )), (21) 

where m# represents the mass of the soliton which is taken as the mass of nucleon, for the 
ground state AQ baryons. 

The binding potential between the heavy rneson and the chiral soliton is simple 
harmonic [22] and hence the wave function is taken as 

(22) 

K is the spring constant and its value is taken to be (440 MeV) 3 [23]. In the rest frame of 
the initial state, v^ = (1,0) and \' 2 directed along z-axis we obtain the Isgur-Wise 
function (21) using (22) for non-relativistic recoils i.e., \\' 2 \ 2 w 2(uj2 1), as 





2 .. ,- (23) 

Again cui is obtained by considering the kinematics of the system, given as 



Now with eq. (24) the baryonic Isgur-Wise function (23) is found to be 

77(^3) =0.263. (25) 

4. Results and discussion 

In order to obtain the value of the ratio of the branching fractions Br(Afc >AJ/*0)/ 
Br(B >//^) with eq. (13) we use the following values. The quark masses are taken as 
m b = 4.5 GeV and m s 550 MeV. The masses of A. h and fi particles are taken from ref. 
[1] as M Afc = 5621 MeV and M B = 5281 MeV. The other particle masses and the life times 
of Ab and of 5 are taken from ref. [24]. With these values we obtain the ratio of the 
branching fractions to be 

/A = 1.017. (26) 

Using the value <r A ,/cr B o =0.1/0.375 from ref. [1], we obtain the production cross 
section times branching fraction for the decay A b -+AJ/ij> relative to that of the decay 
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a271 ' (27) 

which is in excellent agreement with the CDF data 0.27 0.12 0.05 [1]. The overall 
agreement of our result with the experimental data justifies the fact that factorization 
approximation works well for the description of color suppressed nonleptonic decays 
(decays proportional to a 2 ) of the heavy hadrons. 

In this paper we have estimated the ratio of the branching fractions of the color 
suppressed nonleptonic decays A fc ->AJ/^ and fi-> //?/>, based on the factorization 
approximation. Treating the s-quark as heavy, we have used HQET for the evaluation of 
the hadronic matrix elements. Since HQET does not predict the shape of the IW function 
and of course visualizing the fact that the IW functions are different for the baryonic and 
mesonic sectors, due to the different configuration in their light degrees of freedom, we 
employ the well established models for their evaluation. The Isgur-Wise function is 
evaluated in the quark model [8] for B-+KJ/i/j decay and in the bound state soliton 
picture [15] for A*-* AJ/i/> decay mode. Since the evaluation of the Isgur-Wise function 
is of course model dependant, there are some theoretical uncertainties introduced as two 
different models are followed for the estimation of mesonic and baryonic IW functions. 
The ISGW quark model for the mesonic case is quite well established [6-8, 17]. In fact if 
the baryonic IW function could have been evaluated in the quark model then hopefully 
the uncertainties would have been small. However in the absence of such a treatment we 
have followed a different model for it, so it is necessary to test the reliability of the model. 
For the sake of comparison we therefore calculate the branching ratio for B*KJ/ijj 
process using eqs (8) and (20) and thus from eq. (26) we obtain the branching ratio for 
A/, ^AJ/-0 decay process. Now comparing this predicted value with the quark model 
calculation [25] (without using HQET) the validity of the bound state soliton picture for 
the baryonic IW function can be seen easily. 

To calculate the branching fraction for B Q -+KJ/i/j, we use |V CJ | =0.9738, |V cfc | = 
0.038, 2 = 0.23. The decay constant/^, can be obtained from the leptonic width of the 
charmonium state [11] as/ v , = 382 MeV Thus we obtain Br(B->//^) to be 



Br(BP-+KJ/il>) = 3.917 x 10~ 4 , (28) 

and with eq. (26) we obtain 

Br(Az,->AJ/-0) = 3.984 x 10~ 4 . (29) 

Using 1 //WQ corrections to the baryonic form factors in the quark model Cheng et d [25] 
obtained 

Br(A fc -* A//7/;) = 2.1 x 10~ 4 . (30) 

Our present calculation is somewhat greater than their predictions. This decay mode is 
now observed experimentally by CDF [1]. The branching ratio of Br(A* - A//VO = 
(3.7 1.1 0.4) x 10- 4 assuming Br(5 -> KJ/ifi = 3.7 x 10~ 4 . It is seen from eq. (30) 
that the quark model [25] prediction does not agree well with the experiment. In contrast, 
our calculation is very simple and free from any theoretical uncertainties besides in the 
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sector of Isgur-Wise functions. It is interesting to note that our theoretical predictions 
(28) and (29) agree reasonably well with these experimental values. Therefore it should 
be inferred from the above observation that the evaluation of the IW function using two 
different models do not introduce much theoretical uncertainties. Our results indicate that 
the heavy quark symmetry does work in the first approximation to calculate the exclusive 
decays involving b > s transitions. 
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AL. This leads to bounds on Higgs mass. This method has been used by several authors 
[2,6]. 

An alternate and in our opinion more transparent approach to the 'triviality' problem in 
< 4 -theories has been provided by Stevenson using variational methods for calculating the 
effective potential [3]. In Stevenson's approach, the relation between the bare and 
renormalized coupling constant AREN (defined as the fourth derivative of the effective 
potential at the origin and thus different from renormalized coupling mentioned above) is 

[l-moMma)] . } 

AREN = AQ , ... -, rr- , (** 

[1 + 6A /_i (m R )] 

where m| is the second derivative of the effective potential at the origin and 

(5) 



Stevenson showed that in the limit of infinite cut-off, the ground state is unstable (i.e. the 
energy is unbounded from below) for any positive finite AQ. This result goes against our 
naive expectation since the tree level effective potential is bounded from below. However, 
eq. (4) clearly shows that for infinite cut-off AREN = 2Ao, and this is the reflection of the 
energy of the ground state being unbounded from below (see eq. (4.24) of ref. [3]). Then 
for a stable ground state, AQ and AREN both tend to zero for infinite cut-off, which is just a 
trivial theory; we thus arrive at the conclusion identical to the one referred to earlier. 
Once again, if we regard the </> 4 -theory as an effective one, AQ can have value up to some 
maximum below which a stable ground state is possible. For a given coupling constant 
the Higgs masses can be calculated and for the effective theory to be physically 
meaningful we should have all these masses sufficiently below cut-off A c , say A c /5 [5]. 
With such a restriction, one can establish that the Higgs quartic coupling cannot be 
arbitrarily large, i.e., it has a maximum allowed value, which translates into an upper 
bound for MH/MW- 

Triviality bounds in scalar theories have been studied by various non-perturbative 
methods. These include RGE [5], improved perturbative [6], and lattice theoretic approach 
[7]. Results are similar, namely that MH cannot be heavier than a value in the range 
800GeV-lTeV [4]. Recently we have used a variational approach [8] to arrive at a 
similar result. This last approach is extremely simple and admits easy generalization to 
situations more complicated than SM. 

Supersymmetric (SUSY) generalization of the SM have been studied in recent times 
[9]. The most economic SUSY-extension of the SM is a minimal (MSSM) one [9]. In this 
version, the quartic couplings are restricted by the gauge coupling with the result that the 
Higgs cannot be arbitrarily heavy. At tree level, one has the relationship 

M^<M|cos 2 2/3<M^, (6) 

where tan/5 is the ratio of the vacuum expectation values (VEV's) of the neutral 
components of the scalar fields HI and H\ that the MSSM involves. Going beyond the 
tree approximation does not change (6) qualitatively. Thus, Quiros [10] gives the bound 

M H <125GeV (7) 



An alternative supersymmetric model proposed is the next to minimal supersymmetric 
model (NMSSM) which has two SU(2) <g> 17(1) Higgs doublets and one Higgs singlet [9]. 
The inclusion of a Higgs singlet is suggested in many superstring models and grand 
unified models. The NMSSM has more coupling parameters than the MSSM and hence it 
is an interesting theoretical question to enquire into the upper bounds on the Higgs 
spectrum, particularly the lightest of them. We expect this to be much higher than the one 
given in (7) and this is the subject matter of our investigation. 

The method we follow here is a variational one. Starting with Hamiltonian of the 
NMSSM we use a gaussian trial wave functional for the ground state and obtain estimates 
of mass spectra in terms of the bare parameters of the theory. The strategy then is to vary 
bare parameters over their entire range, impose restrictions that the masses cannot get 
very close to the cut-off (say less than Ac/5) and obtain the highest mass of the lightest 
Higgs particle. The parameter space is however very large, and we will be making 
specific choices of parameters in the hope that our results will be typical of the model 
itself. 

In the context of NMSSM, recently Wu [1 1] has analysed the problem (triviality bound 
of Higgs mass) on the basis of one-loop /3 function following the work of Dashen and 
Neuberger [2], imposing the restriction that the heaviest Higgs mass in the theory be 
smaller than the Landau pole position AL. However his estimate of the Landau pole is 
based on one-loop result of the (3 function which is questionable since the coupling 
constant becomes very large in the neighbourhood of the Landau pole. Further the scalar 
masses have all been calculated on the basis of 'tree approximation' which again 
becomes uncertain for large coupling. Our variational method avoids both these 
problems, since the concept of triviality that we use following Stevenson [3] does not 
directly depend on the position of Landau pole. Further we determine all masses by 
variational rather than perturbative methods, which a priori does not restrict the coupling 
to be small. 



2. The model 

In the NMSSM, the potential of the Higgs sector [1 1] is 
V = 



h.c.) +m$N*N + n%(N 2 + JV* 2 ). (8) 



Here ^ = (</>J, </>?) and $ 2 = (<4 $) are two 517(2) <8> 17(1) doublets, N a complex 
singlet, ra's are mass parameters, g's are gauge couplings and h, X are Higgs couplings. 
The last five terms represent SUSY-breaking. Equation (8) has two coupling constants h 
and A; we will study the strong coupling behaviour when h is very large and hence for 
simplicity we set A = 0. Also we take m 3 to be real and m, = m 2 = m for simplicity. 
It is more convenient to work with the fields defined by 
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Now the Higgs potential reduces to 

V = h 2 \N\ 2 ( X \X\ + X 2 X2) + \h 2 \x\Xi - X 2 X2 + xlxi ' X\X2\ 2 

+ (m 2 - m 2 )x{ Xl + (m 2 + m 2 } X \X2 + n$N*N + m 2 (N 2 + AT 2 ). (10) 

We now assume that the fields xi,2. where 



(the superscripts c and denoting charged and neutral components), break the 517(2) 
U(\] symmetry by assuming a non-zero vacuum expectation value (VEV) 



All other fields in (10) are assumed to have zero VEV's. Writing N = \/\/2(N\ + iN 2 ], 
there are ten real fields in (10): XIR, Xw X Q \ R > Xw X 2R , X 2[ , X 2R , X/> M and ^2, which we 
denote by 77,- (i = 1 to 10) respectively and their tree level masses by M,-. Even after xi 
develops a non-zero VEV, the Higgs potential has a residual symmetry. To see this 
explicitly we define: 

(%, -m,74), 



5 i = X = % and S 2 = Xy = ^8- 
Now shifting the fields by their VEV's, Higgs potential can be written as 

V =\h 2 (N 2 +N 2 )(G 2 + (S } + v) 2 + H 2 + S 2 ] 

+ i/z 2 [G 2 + (Si + </) 2 - H 2 - 5 2 ] 2 + \h 2 [G H - (5! + t;)5 2 ] 2 
+ i (' 2 - !) [G 2 + (5i + u) 2 ] + i (m 2 + m 2 ) (H 2 + S 2 } 
+ 1 (m 2 + 2m 2 )N 2 + 1 (m 2 - 2m 2 )N 2 . 

Under a 517(2) rotation wherein G and H are triplets, and (Si ,S 2 ,Ni, N 2 ) all singlets, V is 
invariant. We expect then M\ = M 2 = M$ and MS = Me = M-j on account of this 
symmetry. We note that this symmetry is present in a two-doublet model for the special 
choice of our parameters. Thus in the two-doublet model discussed in 2 of ref. [12], 
setting v-i = v 2 = = and AI = A4 makes H^, H^ and H degenerate. We will build in 
this residual symmetry in our variational approach. 
In terms of ten real fields the Higgs potential 

V =(m 2 - m 2 ), 2 +ifcV +i(m 2 - m| +itfv%? + tg +r) 
+ i (m 2 - m 2 + -h 2 v 2 }r, 2 + - (m 2 + m 2 - ^h 2 v 2 } (r/ 2 + r, 2 + 77?) 
+ -(m 2 + m 2 -I- -h 2 v 2 }r,l + - (m 2 + 2m 2 + /z 



7^1 cn vr n ^ i?ak...r mno 



Higgs mass in NMSSM 
+ l(mJ-2 

+ o h 2 [ril(r) 
o 



7% - rft 



i, + il ~ V - ~ T) + v(n ~ r - r} - r + r, 2 ) 



+ r]l) + 7$$] + linear + Vcubic, (13) 

where V\\ neai and Vc U bi c respectively represent the terms linear and cubic in fields. The tree 
level minima condition is equivalent to equating the term linear in 773 (i.e. X\R) in (13) to 
zero, 

m 2 ~m\ +i/zV = 0. (14) 

Inspection of the quadratic terms in (13) together with (14) immediately tells us that 

M 1 =M 2 =M 4 = 0, (15) 

indicating that 771, 772 and 774 are the Goldstones. Furthermore, we have for the charged 
Higgs triplet 



i/rV, (16) 

where the degeneracy is as expected. Lastly, 

M| =ra 2 -m^+|/iV, (17) 

M 2 =m 2 + m 2 + i/zV, (18) 

Ml = m\ + 2m 2 + h 2 v 2 (19a) 

and 

M 2 = m 2 4 - 2m] + h V (19b) 

are the masses of neutral Higgs. For simplicity we also take m*, = so that 

Ml = M 2 . (20) 

3. Gaussian trial wave-functional 

In order to obtain information on masses beyond the tree level, we follow a variational 
method with a Gaussian trial wave-functional. Most generally, this wave functional 
would be the vacuum state of a set of free fields of masses Oi , . . . , ftio with the Q's (and 
v in eq. (12)) representing the variational parameters. However, taking a variational 
ground state that respects the residual symmetry stated in the last section, we set 
Qj = fi 2 = ^4 and f^ = Og = ^7- Further imposing the symmetry in (20) we also put 
Op = IO = HAT. We thus have five independent masses fii, f^, ^5, Qg and fi#, and of 
course v, as variational parameters. 
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Following standard techniques [13], the expectation value of Hamiltonian density H. in 
our trial vacuum state wave functional is 



V G =(m 2 - 



-h 2 v 4 + 2 (M 



Here 



and 



+ 37 (Q 2 )(7 (Q 2 ) - 

= L_ / d 4 fe ln(fcf + Q 2 ) + constant 



Differentiating (21) with respect to fijy, fi 2 , 



= M 
= M 



(22a) 



fi 2 - M 2 + |/i 2 [4/o(n 2 ) - 3/ (0 2 ) 
Using equations (23), VG reduces to 



(22b) 

and fig gives us five mass equations: 

(23a) 
(23b) 

, (23c) 
(23d) 

. (23e) 



+ [27, (O 2 ,) + 37! (fi 2 ) + h (Of) + 37: (ft 2 ) + 7! (H 2 
- i^ 2 [57 2 (n 2 ) .+ 7 2 (0 2 )+ 57 2 (H 2 ) + 7 2 (0 2 )] 



2 )7 (^)]. (24) 

Here O/'s are to be understood as depending on v through (23). Differentiating V G with 
respect to v 2 , we get 



2 



(25) 



Higgs mass in NMSSM 

Setting dVa/dv 2 to zero and using (23c), we get 

Sll = h 2 v 2 . (26) 

It is clear from (25), that increasing h would increase Sl\ without limits. However, the 
limit would be set by demanding that stationary solution (25) be stable, i.e. the stability 
matrix (<9 2 Vo) be positive definite. Stability condition is obtained by considering (d 2 Vo) 
to be a function of six variables five mass parameters (O 2 , ft 2 , fi 2 , fig and fi 2 ^) and v 2 , 
and demanding all its eigenvalues to be positive. As in the case of standard model [8], we 
expect that as h increases, this condition would no longer be satisfied beyond a certain 
maximum value of h. 

Our query regarding triviality bounds does not involve the complete numerical 
solution. Of five independent masses ^1,^3,05,^8 and 1 N , O] is the mass of the 
Goldstone bosons. This in an exact calculation is expected to be zero but in variational 
methods (see ref. [13] for elaboration), we can get a small but non-vanishing mass. Of the 
remaining four Higgs masses, we wish to find out whether they can be made arbitrarily 
heavy relative to v 2 (or A/^). The only condition we would impose is the same as laid 
down by Hasenfratz and Nager, namely that for a cut-off theory to make any physical 
sense, each one of the masses fZ/ must not be close to or greater than the cut-off ; we put 
an upper limit of A c /5 for defmiteness, for all fi/'s (see ref. [5]). 

Our task is then to set $Vs at me i r maximum possible values and determine the value of 
h for which the eigenvalues of stability matrix go from positive to negative. Since there 
are only three mass input parameters m, m 3 and 1114, all the five fi/'s cannot be assigned 
arbitrary values by suitably choosing m's. Furthermore, only O^ involves m 2 , so that we 
can set fl w as the highest acceptable mass namely A c /5 right away. Of the remaining four 
masses, we immediately have the sum rule, 

a Q u Lc - QUO Jui 

2 ) + 7 (f2 2 ) - I (tf s )]. (27) 



Also since /o(0 2 ) is a decreasing function of O 2 , we get from (23) and (27) 

til > ft 2 and fi 2 , > n 2 . 
The first possibility is to set 

f} 8 = f2 3 = Ac/5; 

QI is the Goldstone and 1$ is then the lightest Higgs mass. However in this case f2 5 = fii, 
and the maximum value of O 5 turns out to be lower than in the next sequence of masses. 
Next we assume Q 8 = A c /5 and the sequence 



In this f2 3 is the lightest Higgs. From (27), pushing ^5 towards fig will make n 3 go 
towards HI, i.e. will make Q 3 lighter. We then expect some kind of optimal situation to 
arise if fi 5 = 3 . With the values of 3 (obtained from (23)) corresponding for various h, 



Using this value, then the upper bound on mass of the lightest Higgs MLR = f^ = LZs is 

_ in 1 
IU.J. 



The Goldstone mass l\ for this choice, as we stated before, is not zero but smaller than 
other masses in the spectrum. 

4. Conclusion 

We have shown that in a certain range of parameters of the NMSSM, the triviality upper 
bound on the mass of the lightest Higgs particle is ~ lOMw- Relaxing restrictions on the 
parameters can lead to even larger masses for the lightest Higgs. We have however not 
attempted to determine the absolute bound taking into account the full range of variation 
of parameters in the NMSSM (including SUSY breaking parameters). This is because our 
main aim was to show that in the non-minimal version of supersymmetric model, one 
need not be constrained by the rather strict limits on the Higgs mass that one obtains in 
the minimal model. 

We would also like to comment on the upper bound of about 140 GeV for the lightest 
Higgs boson mass obtained by several authors [14]. This arises by the requirement that 
the SUSY theory remain perturbative up to some scale MQUT ~10 15 GeV. This clearly 
forces the quartic coupling constant to remain small. From a phenomenological point of 
view where the couplings of the NMSSM model have to be realistic over a huge energy 
range, such restrictions are eminently reasonable and hence in that respect the bounds 
obtained by reference [14] are phenomenologically more relevant. We, in our analysis, 
have not constrained the coupling constant from any such requirement and thus the limits 
obtained are much higher. Comparing our result with that of Wu [11] who obtained a 
much lower bound by a factor of four, we find that the criterion used by him to define 
triviality bound (as explained in the introduction) is different from ours. Also the method 
of calculating the scalar mass particles is based on tree approximation unlike our non- 
perturbative variational method which goes beyond the tree approximation. The difference 
in our numerical bounds we believe can be attributed to these differences in our 
approaches. Also limits on the Higgs mass above 1 TeV are of little interest. There is no 
possibility in any near future to detect any signals for such a heavy Higgs. Moreover 
Higgs particle with masses above 1 TeV with widths comparable to masses will make the 
mass parameter rather meaningless from an experimental point of view. Further, it is 
interesting to note that the bound obtained for the lightest Higgs (~ 10 M w ) is of the same 
order as the bounds in the SM suggesting that in a supersymmetric theory, the Higgs mass 
bounds have the same features as the regular SM. 
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Abstract. Excitation functions and a few isomeric cross-section ratios for production of (1) 192 Au, 
193 Au, 194 Au, 195 Au and )92 Ir nuclides in a-induced reactions on 191 - I93 Ir, (2) ""TI, 197w Hg, 198 '>*T1, 
199 T1 and 200 T1 nuclides in a-induced reaction in 197 Au and (3) 183 Re and 184m 'Re nuclides in as- 
induced reaction in 181 Ta and 185 Re are obtained from the measurements of the residual activities by 
the conventional stacked-foils technique from threshold to SOMeV. The excitation function and 
isomeric cross-section ratios for nuclear reaction 181 Ta(a, n) 1 "'"^Re are compared with the theore- 
tical calculation using the code Stapre which is based on exciton model for pre-equilibrium phase and 
Hauser-Feshbach formalism taking angular momentum and parity into account for the equilibrium 
phase of the nuclear reaction. All other experimental excitation functions are compared with the 
calculations considering equilibrium as well as pre-equilibrium reaction mechanism according to the 
geometry dependent hybrid (GDH) model and hybrid model of Blann using the code Alice/91. The 
high energy part of the excitation functions are dominated by pre-equilibrium reaction mechanism 
whereas the low energy parts are dominated by equilibrium evaporation with its characteristic peak. 
The GDH model provides a potentially better description of the physical process (i.e., a higher pro- 
bability for peripheral collisions to undergo precompound decay than for central collisions) compared 
to hybrid model. However in the energy range of present measurement most of the excitation func- 
tions are fitted reasonably well by both GDH model and hybrid model with initial exciton number 
N Q = 4(N n = 2, N p = 2, N h 0). Barring a few reactions we have found the overall agreement 
between theory and experiment is reasonably good taking the limitations of the theory into account. 

Keywords, a-induced reactions on Ir, Au, Re and Ta; stacked foil technique; isotope production; 
equilibrium and pre-equilibrium decays; Stapre and Alice/91 codes. 

PACSNo. 25.60 



1. Introduction 

It is well known that the compound statistical model coupled with the exciton model gives 
a correct overall description of the excitation functions and particle energy spectra in 
nuclear reactions at medium energies. However, the calculations partially fail to account 
for details such as the exact position of the maximum or the slope of the ascending and 
descending part of the excitation functions. The high energy parts of the excitation 
functions are dominated by pre-equilibrium reaction mechanism whereas the low energy 
parts are dominated by evaporation with its characteristic peak. The pre-equilibrium 
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processes have been investigated by several authors [1-6]. Excitation functions for 
equilibrium and pre-equilibrium reactions have been studied in (a, xn) and (a, pxn) 
reactions [7-10]. Theoretical calculations for pre-equilibrium process have been carried 
out in terms of exciton models [1-3]. Several models [1-6, 11-15] have been proposed to 
explain the emission of energetic light particles by the equilibration process from the 
nuclear system excited at medium energies. Predictions from these models with regard to 
excitation functions and the energy spectra of the emitted particles compared well with the 
existing experimental data. This has prompted a continued interest in these models as tools 
both to predict cross-sections for a number of nuclear reactions for practical purposes and 
to test the adequacy of the underlying physics. The isomeric cross-section ratios for a pair 
of isomeric states are known to depend strongly on the spins of the isomers concerned as 
well as on the spins of the higher lying levels populating the isomers. Experimental and 
theoretical studies on the isomeric cross-section ratios specially as function of incident 
particle energy should therefore lead to useful information on the spin cut off parameter as 
well as on the level structure of the residual nuclei. The excitation functions and isomeric 
cross-section studies are significant for the investigation of the mechanism of nuclear 
reactions. The present work on alpha induced reactions on the target nuclei 181 Ta, 185 Re, 
191 Ir and 193 Ir is also intended to supply some new data in the alpha energy range from 
threshold to 50.0 Me V. Since the nuclear reactions 197 Au(o;, xnypza) have already been 
studied extensively by several authors (see [16] and references therein), we have measured 
the cross-sections and the isomeric cross-section ratios for 197 Au(o:, xnypza) reactions so 
as to standarize our measurements by comparing our measurements with previous 
measurements. Hence the nuclear reactions 197 Au(a, xnypza) have not been given much 
importance for theoretical analysis. However the nuclear reaction cross-sections are being 
presented in this paper only for completeness. 

This paper continues a series of reports on gross features of the interactions of 
intermediate energy light projectile with medium and heavy nuclei. The experiments 
have been performed at the Variable Energy Cyclotron Centre, Calcutta. Several 
excitation functions for the reactions 191 ' 193 Ir(o!, xn) I95 ~ x > 197 ~ j: Au, 191 Ir(a, an) 192 Ir, 
197 Au(a, xn) 20 -*T1, 197 Au(a, p3n) 197 Hg, 185 Re(a, an) 184 Re and 181 Ta(a, xn) 185 -*Re 
measured by the stacked foil technique are presented. The excitation functions of the 
radioactive products observed in reactions contain some information about the 
mechanism of the interaction of the a-particles with Ir, Ta Re and Au nuclei. In this 
work calculations in the framework of the equilibrium statistical model and pre- 
equilibrium model using the code Alice/91 [17] and Stapre [18] were performed and the 
results are compared with experimental excitation functions. The Stapre code [18] which 
is based on exciton model for pre-equilibrium phase and Hauser-Feshbach formalism 
taking angular momentum and parity into account for the equilibrium phase has been 
used to calculate isomeric cross-section ratios and the results are compared with the 
measured isomeric cross-section ratios. 



2. Experimental procedure 

Excitation functions for alpha induced reactions on 181 Ta, 191 Ir, 193 Ir and 197 Au were 
determined using the absolute yields of characteristic 7-rays pertaining to the decay of 
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each radioactive residual nuclide as usually done in stacked foil technique. The targets 
were commercially available thin self-supporting foils. The thicknesses of gold foils were 
4.837 |im, that of irridium foils 13.021 urn, and tantalum foils were 10.521 urn obtained 
by weighing. Then the stacks were irradiated in a chamber (as shown in figure 1 of [9]) 
specially constructed for this purpose having the facility to irradiate up to four stacks one 
after another. The beam spot on the targets was limited to 5.0mm in diameter by using a 
10.0 cm long aluminum collimator in front of the targets. The Au and Ta stack consisted 
of about 20 targets whereas Ir stack consisted of about 15 targets. Since the Au foils were 
very thin, these were interspersed with 23.4 urn thick Al degrader-cum-catcher foils. The 
stacks were exposed to the unanalysed external beam from the 224 cm variable energy 
cyclotron in Calcutta. The beam current on the targets was kept below 200 nA. The total 
<a-particle beam was collected and measured using a calibrated ORTEC current 
integrator. The absolute cross-sections for the reactions 197 Au(a, xnypaz) were also 
compared with previous measurements [16]. The cross-section for the a-induced reaction 
on Au are very well known [16]. Comparison of 200 T1, 199 T1 and 198 T1 production cross- 
section from previous measurements with our measurements thus provides high 
reliability to the absolute cross-section measurements. The unanalysed beam energy 
resolution was ~ 0.2 MeV. The accuracy in absolute energy is expected to be = 2.0 Me V. 
However from some of the excitation functions for the a-induced reactions on Au the 
energies were re-calibrated with respect to their thresholds and thereby improving the 
accuracy in absolute energy to = 0.5 MeV. 

The mean beam energy at the half-thickness of each target of a stacked foil assembly was 
calculated from the energy degradation of the initial beam energy according to the given 
stopping power values for different materials. We have used the stopping powers from the 
tabulated values of Williamson et al [19]. However the tabulated values are in multiple of 
2.0 MeV up to 50.0 MeV and in steps of 5.0 MeV thereafter. Therefore interpolation is 
required at every foil of the stack. To improve the accuracy of interpolation the stopping 
powers S(E] are fitted by a non-linear least square method to a function of the form: 



S(E] = A, E>- ] +, B 

1=1,4 =1,2 

where A,-, Bk and Q are the constants to be determined from non-linear least square fit. 
The fit improves the overall accuracy of the stopping power interpolation. The average 
thicknesses of the target and degrader foils were determined by weighing. Each foil was 
cut out into a square shape and pasted over an annular Al ring having 30.0mm as the 
outer diameter. 

The 7-rays emitted by the activated foils were detected either with a Ge(Li) detector or 
with a HPGe detector (both having 30% efficiency) available at our centre. In most of the 
cases the 7-rays, used in yield determination as listed in table 1, stand out very 
prominently in the spectra and does not pose any identification problem. The 7-ray 
spectra from the Ge(Li) and HPGe detectors were analysed using personal computer 
based multi channel analyser and were stored on hard disk or floppy diskettes. The Super- 
32 computer was used to analyse the 7-ray spectra stored on these hard disk and floppy 
diskettes. One program EXPANL1 was written to select the required peaks from the 7-ray 
spectra stored on hard disk and floppy diskettes and prepare a data file EXP2DT. Another 
program ISABELLE written by Kern [20, 21] (modified and renamed as EXPANL2) was 
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Table 1. Half-lives, 7-energies, branching ratios of the 7-decays and Q-values for a- 
induced reactions on Ir, Au, Re and Ta. 



Juclide 


Half-lives 


E 7 (keV) 


!,(%) 


Reaction 


Q-values (MeV) 


n Au 


4.94 h 


316.51 


58.00 


191 Ir(a : , 3n) 


-25.709 










193 Ir(a,5n) 


-39.679 


33 Au 


17.65 h 


186.17 


10.11 


I91 Ir(a,2n) 
193 Ir(a,4n) 


-16.944 
-30.914 


94 Au 


38.02 h 


328.50 


60.00 


191 Ir(ct,n) 


-10.084 










193 Ir(a,3n) 


-24.054 


95 Au 


1 86.09 d 


98.88 


10.90 


193 Ir(a, 2n) 


-15.670 


92 Ir 


73.83d 


316.51 


82.81 


I93 lr(a,2p3) 


-36.053 










!93 Ir(a an) 


-7.757 










19I Ir(a,2pn) 


-22.084 


97 T1 


2.84h 


425.84 
152.23 


12.91 
07.21 


197 Au(a,4n) 
197 Au(a,4n) 


-32.606 
-32.606 


97 "'Hg 


23.80h 


133.96 


34.00 


197 Au(a,p3n) 


-29.678 


98g T1 


5.30 h 


411.80 


82.10 


I97 Au(a,3) 


-25.434 


9Sm T j 


1.87h 


282.80 


28.00 


197 Au(a,3n) 


-25.434 


99 T1 


7.42h 


455.46 


12.31 


197 Au(a,2) 


-16.813 






208.21 


12.21 


I97 Au(a,2n) 


-16.813 


:oo T1 


26.10 h 


367.94 


87.20 


197 Au(a,n) 


-09.736 


83 Re 


71.20d 


162.33 


23.58 


18t Ta(a,2n) 


-16.340 


84 Re 


38.00 d 


792.07 


37.64 


181 Ta(a,) 


-09.873 


84m Re 


169.8 d 


252.85 


10.89 


181 Ta(a,n) 


-09.873 


84 *Re 


38.00 d 


792.07 


37.64 


185 Re(a,2p3n) 


-35.975 


l8 4,n Re 


169.8 d 


252.85 


10.89 


185 Re(a, 2p3n) 
I85 Re(a,an) 


-35.975 
-07.679 



h = hour, d = days, y = year. 

used to analyse the peaks stored in the data file EXP2DT. The program ISABELLE fits 
the 7-ray spectra very well (see figure 2d in [9]). The integration of the analytical 
function as well as experimental peak was used for yield determination. The 7-ray yield 
was determined by integrating the experimental peak points and subtracting the 
calculated background from it. 

The efficiency calibrations of the detector were done with a standard 152 Eu and 133 Ba 
radioactive source obtained from B.R.I.T. (Department of Atomic Energy) available at 
our centre. The efficiency of the detector was interpolated to the required energy value 
from the measured efficiency curve. However to improve the interpolation, the efficiency 
curve was similarly fitted by a non-linear least square method to a function of the form 
given in (1). The range of the fit was from 122keV to 1408keV. The fit improves the 
accuracy of interpolation considerably. For 7-ray having energy 7 < 122keV, the 
efficiency of the detector was extracted by manually plotted efficiency curve using low 
energy 7-rays from 152 Eu and I33 Ba calibrated sources. 

The nuclear data necessary for the evaluation of the cross-sections are presented in 
table 1. The half-lives of the radioactive atoms are taken from the chart of nuclides [22], 



the 7-ray energies and branching ratios are taken from the table of isotopes [22]. In 
table 1, only those 7-ray s are listed which were chosen for the calculation of the cross- 
sections. Also included in table 1 are reaction ^-values. <2-values were calculated using 
the atomic mass table of Wapstra and Audi [23]. 



2.1 Cross-section determination 

The number of observed decays Z per unit time is related to the total number of decays ZQ 
per unit time by 



Z = Z/[e(E 7 ) * 7 7 (abs)], 



(2) 



(a) 



193-IR(ct.2n)195-Au 
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Z 
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(b) 



193-Ir(a, 3n)194-Au 
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(d) 
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Figure la-d. The total residual production cross-section in mb for the reactions (a) 

193 Ir(a, 2) 195 Au() (b) 193 Ir(<*, 3rc) 194 Au (),((:) t 193 )' 191 !^, (4)2] 193 Au[( B ),a] 
and (d) 191 Ir(o!, 3n) Au() are plotted as a function of a-particle bombarding 
energy. The solid line is the geometry dependent hybrid model and the circle-line (-Q-) 
is the hybrid model fits using the code Alice/91 with NQ = 4,JV n = 2,N P 2 and 



7i en vr. 



^.,^ M , tone 
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where e( 7 ) is the detector efficiency and 7 7 (abs) is the absolute 7-ray abundance per 
decay. The total number of decays, ZQ, is related to total reaction yield R for simple 
decays (simple decays correspond to direct production of radio isotopes by the nuclear 
reactions and we have used only such decays in all the measurements reported here) by 
the relation 

R = Z exp(Ar 2 )/[(l - exp(-At)) (1 - exp(-Ar 1 ))/Ar 1 ], (3) 

where A is the decay constant, T\ the irradiation time, TI the cooling time after the end of 
irradiation and the beginning of the measurement, T the counting time and R is related to 
the cross-section a by the relation 

5 X }-J], (4) 



Table 2. Experimental cross-section for oi-induced reaction with 191 Ir 
and 193 Ir. 


Energy in 
(MeV) 




From 19] Ir 




192-Au(3n) 


193-Au(2n) 


194-Au() 


48.4 
0.9 


580.62 
0.30% 






45.8 
1.0 


656.41 

0.70% 






43.8 
1.0 


727.27 
0.80% 






41.9 ' 
1.0 


939.84 
0.70% 






39.7 
1.1 


1141.51 
0.40% 






37.6 
1.1 


1562.89 
0.30% 






35.4 
1.1 


1602.75 
0.30% 


110.9* 
0.80% 




33.1 

1.2 


1349.15 
0.40% 


236.26 
1.50% 




30.7 

1.2 


841.48 
0.60% 


438.09 
2.50% 


23.3* 
0.15% 


28.1 

1.3 


172.789 
2.10% 


640.60 
1.40% 


29.8* 
0.24% 


25.4 
1.4 




522.11 
0.30% 


32.7* 
0.80% 


22.6 

:.f 1 .5 




269.65 
0.30% 


35.6 
1.60% 


10.5 
1.6 




45.74 
0.70% 


33.64 
1.00% 


16.1 
1.8 






3.14 
1.00% 
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Table 2. (Continued) 



Energy in 
(MeV) 


Cross-section of the product in 


milli bams from 193 Ir 




195-Au(2n) 


192-Ir(an) 


194-Au(3n) 


193-Au(4n) 


48.4 
0.9 


47.06 
8.20% 


196.998 

0.46% 


179.462 
1.30% 


854.336 
2.00% 


45.8 
1.0 


49.72 

5.27% 


83.863 
0.46% 


207.039 
1.10% 


899.380 
3.00% 


43.8 
1.0 


57.23 
5.30% 


52.057 
1.00% 


320.463 
0.90% 


1077.775 
1.00% 


41.9 
1.0 


64.28 
4.00% 


37.069 
0.69% 


445.464 
0.60% 


904.602 
0.80% 


39.7 
1.1 


59.18 
5.00% 


28.935 
1.50% 


608.938 
0.30% 


644.735 
0.60% 


37.6 

1.1 


84.80 
2.60% 


25.356 
1.70% 


1028.389 
0.16% 


441.116 
0.60% 


35.4 
1.1 


110.63 

2.54% 


18.491 

1.75% 


1342.906 
0.08% 


118.3* 
1.80% 


33.1 

1.2 


139.56 

2.07% 


11.801 

3.07% 


1318.248 
0.12% 




30.7 
1.2 


303.44 
0.96% 


6.505 

6.10% 


1061.1 
0.15% 




28.1 
1.3 


580.44 
0.52% 


4.193 
9.10% 


478.3* 
0.24% 




25.4 
1.4 


757.07 
0.42% 


2.053 
9.90% 


17.0* 
0.80% 




22.6 
1.5 


426.80 
0.93% 


1.686 

10.0% 






19.5 
1.6 


97.49 

2.73% 


1.686 
10.0% 







where N A is the number of atoms/cm 3 of the target material, 8x is the thickness of the foil 
(in cm) and J is the total number of incident particles during the irradiation (calculated 
from total charge measured by current integrator). 



3. Experimental results 

In table 2 and figures (la) to (Id) our experimental results for the production of 192 Au, 
192 Ir, 193 Au, 194 Au and 195 Au radio-nuclides via ^-induced reactions on 191 . 193 Ir are 
summarized and plotted as a function of a-particle bombarding energy. The experimental 
cross-sections are in millibarns and presented along with the percentage of statistical 
errors only. Since the radio-nuclides 193 Au could be produced by I91 Ir(o;,2) and 
193 Ir(a,4) reactions and the radio-nuclides 194 Au could be produced by I91 lr(a,n) and 
193 Ir(o:, 3n) reactions, in the overlapping regions the cross-sections for these reactions are 
taken as proportional to their theoretical values based on the GDH model calculation [6]. 
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Figure 2. (a) The total residual production cross-section in mb for the reaction 
181 Ta(a, n) 184m ' g Re [m=> g, n => w] are plotted as a function of a-particle 
bombarding energy. The solid line is the geometry dependent hybrid model fit using 
the code Alice/91 with N Q = 4,N n = 2 } N P = 1 and N h = 0. The -- and -o- are 
Stapre code calculation with FM = 1000.0 (MeV) 3 and FM = 750 (MeV) 3 respec- 
tively, (b) The isomeric cross-section ratios for the isomer pair 184m 'Re are plotted () 
as a function of a-particle bombarding energy. The solid lined curve is Stapre code 
calculation with FM = 1000 (MeV) 3 . The curve with FM = 750(MeV) 3 is 
indistinguishable with the plotted solid curve. 



(0 



181-Ta(a,2n)183-Re 




30 40 
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Figure 2c. The total residual production cross-section in mb for the reaction 
181 Ta(a, 2n) 183 Re () are plotted as a function of a-particle bombarding energy. The 
solid line is the geometry dependent hybrid model fit using the code Alice/91 with 
N Q =4,N n = 2,N p = 2 and N h = 0. 



These points are marked with symbols (*) in table 2 and beyond these marked points the 
contributions of (a, n) and (a, 2n) reactions to (a, 3n) and (a, 4n) reactions respectively 
are negligible (< 2%). The number of such points for (a, 2n) and (a, 4n) reactions is 
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Table 3. Experimental cross-section for a-induced reaction with 181 Ta. 



Energy in 
(MeV) 


Cross-section of the product in milli barns 


183-Re 


184-Refe) 


184-Re(m) 


*/(*, + *) 


18.8 
1.0 


6.09 

1.5% 


14.45 
0.63% 


11.07 

3.0% 


0.434 


20.7 
1.0 


114.23 
0.3% 


33.50 
0.69% 


28.76 

2.6% 


0.462 


22.6 
0.9 


375.89 
0.09% 


27.23 
0.48% 


23.36 

2.1% 


0.462 


24.2 
0.8 


581.04 
0.11% 


15.78 
1.08% 


14.17 
6.3% 


0.473 


26.2 
0.8 


730.62 
0.08% 


11.62 
0.11% 


9.50 

7.4% 


0.450 


27.8 
0.8 


772.58 
0.05% 


9.07 
0.81% 


8.38 

7.3% 


0.480 


29.4 
0.8 


698.82 
0.11% 


9.13 

1.52% 


6.41 
10.0% 


0.412 


30.9 
0.8 


457.57 
0.13% 


7.54 
1.69% 


7.03 
9.4% 


0.482 


32.5 
0.8 


334.07 
0.15% 


7.65 
1.62% 


6.43 

7.8% 


0.457 


33.9 
0.7 


224.18 
0.18% 


7.14 

1.72% 


7.36 

8.7% 


0.508 


35.3 
0.7- 


177.45 
0.20% 


6.87 

1.73% 


6.73 

7.9% 


0.495 


36.7 
0.7 


132.90 
0.24% 


5.86 

1.91% 


5.87 
5.8% 


0.500 


38.0 
0.7 


118.90 
0.26% 


5.48 
2.03% 


6.11 

7.3% 


0.527 


39.3 
0.7 


91.25 
0.28% 


4.47 
2.20% 


4.97 
9.2% 


0.526 


Table 4. Experimental cross-section for cv-induced reaction with 


185 Re. 






Cross-section of the product in rnilli barns 


(MeV) 


184-Re(g) 


184-Re(m) a, tll 


/fo + ffm) 


32.6 2.0 3.83 5.0% 
36.6 2.0 4.29 4.8% 
40.5 2.0 12.11 1.8% 
44.2 1.8 18.33 0.45% 
47.6 1.7 26.15 0.55% 


4.07 10.0% 
5.28 10.0% 
11. 67 9.0% 
18.73 3.3% 
28. 00 4.2% 


0.515 
0.552 
0.491 
0.505 
0.517 



only one and for (a, n] and (a, 3n) reactions three. Since no data exist in literature for a- 
induced reactions on 191 > 193 lr 5 comparison could not be made with other measurements. In 
table 3 and figures 2a and 2c the total residual production cross-section in mb for 



Table 5. Experimental cross-section for or-induced reaction with 197 Au. 





Cross-section of the product in milli barns 


Energy in 












CMeV) 


200-T1 


199-T1 


198-T10?) 


198-Tl(m) a 


m/K + CTg) 


52.0 0.3 


1.95 1.9% 


37 .29 3. 8% 


174. 15 0.3% 


254.1 8 8. 6% 


0.594 


50.7 0.3 


2.16 1.7% 


45. 12 3.0% 


21 1.19 0.3% 


328.57 7.1% 


0.608 


49.4 0.3 


2.18 1.7% 


44.50 3.1% 


232.51 0.4% 


337.97 5.9% 


0.592 


48.0 0.3 


2.43 1.8% 


47.79 2.3% 


289.24 0.2% 


433.42 3.3% 


0.600 


46.6 0.3 


2.52 1.8% 


5 1.45 1.6% 


364.37 0.2% 


521. 26 2.6% 


0.589 


45.1 0.3 


2.43 1.7% 


48.96 1.4% 


406.23 0.2% 


735.08 2.7% 


0.644 


43.6 0.3 


2.88 1.6% 


57.17 0.8% 


564.90 0.1 8% 


990.33 1.8% 


0.637 


42.1 0.3 


3.08 1.7% 


65.55 1.2% 


750.94 0.3% 


1230.85 1.4% 


0.621 


40.6 0.3 


3.82 0.5% 


82.25 0.8% 


959.55 0.2% 


1417.20 1.2% 


0.596 


38.8 0.3 


4.40 0.6% 


10 1.29 0.9% 


1059.09 0.2% 


1518.60 1.0% 


0.589 


37.2 0.4 


4.90 0.6% 


129.20 0.56% 


101 1.62 0.2% 


1485.62 1.0% 


0.595 


35.5 0.4 


5.50 0.6% 


184.27 0.5% 


925.02 0.2% 


1266.36 0.8% 


0.578 


33.7 0.4 


7.49 0.5% 


317.88 0.3% 


783.06 0.3% 


1023.21 1.3% 


0.567 


31.80.4 


7.03 0.7% 


485. 13 0.2% 


483.44 0.4% 


61 2.74 1.9% 


0.559 


29.9 0.4 


8.28 0.5% 


723.71 0.1 6% 


196.75 0.4% 


208.55 0.6% 


0.515 


27. 8 0.4 


9.32 0.2% 


710.79 0.16% 


23.65 2.2% 


15.44 2.0% 


0.390 


26.5 0.4 


9.61 0.22% 


708.60 0.4% 








25.6 0.4 


10.58 0.4% 


649.1 2 0.4% 








25.6 0.4 


12.53 0.8% 


57 1.78 0.2% 








24.8 0.5 


12. 10 0.2% 


554.83 0.5% 








24.0 0.5 


15. 19 0.4% 


500.91 0.4% 








23.3 0.5 


16.86 0.5% 


345.53 0.33% 








23.1 0.5 


19.90 0.9% 


362.57 1.0% 








22.2 0.5 


22.32 0.2% 


282.92 0.9% 








21.0 0.5 


26.88 0.8% 


173.76 0.86% 








20.8 0.5 


27.84 0.1% 


11 9.30 1.2% 








20.3 0.5 


27.84 0.6% 


90.50 0.8% 








19.4 0.5 


20.82 1.0% 


29.79 0.9% 








18.4 0.5 


10.73 0.5% 


9.41 6.4% 








18.0 0.5 


11. 78 0.1% 


5.57 0.6% 








17.4 0.6 


3.29 0.2% 


0.38 2.4% 








16.4 0.6 


0.54 0.5% 










15.4 0.6 


0.06 7.9% 












Cross section in milli barns 


Energy in 
(MeV) 


197-T1 


197-Hg(m) 








52.0 0.3 


11 60.30 2.0% 


68.85 1.2% 








50.7 0.3 


1242.39 1.7% 


65. 10 1.3% 








49.4 0.3 


1242. 13 2.0% 


6 1.10 2.0% 








48.0 0.3 


1176.41 2.1% 


57. 17 2.0% 








46.6 0.3 


1016.74 2.1% 


54.88 2.2% 








45.1 0.3 


861. 25 2.3% 


5 1.08 2.2% 








43.6 0.3 


629.96 3.8% 


39.11 2.3% 








42.1 0.3 


386.70 6.6% 


25.53 4.4% 








40.6 0.3 


170.06 8.2% 


12.40 5.0% 








38.8 0.3 




2.93 3.0% 








37.2 0.4 




1.81 3.9% 








35.5 0.4 




0.85 6.8% 









the reaction 181 Ta(a, n) 184m '*Re and 181 Ta(a, 2n) 183 Re are summarized and plotted as 
a function of a-particle bombarding energy. The experimental cross-sections are in 
millibarns and presented along with the percentage of statistical errors only. In table 3 and 
figure 2b the experimental results on the isomeric cross-section ratios for the isomer pair 
i84w, SRe by the react i on 181 Ta(or,n) 184m ' s Re are given as a function of a-particle 
bombarding energy. In table 4 the production cross-section as well as the isomeric cross- 
section ratios for the isomer pair 184m 'SR e by the reaction 185 Re(a, an) 184m ' g Re are 
presented as a function of a-particle bombarding energy. Since the isomeric cross-section 
ratios are measured for the first time, no comparison could be made. In table 5 and 
figures 3a-e the total residual production cross-section in mb for the reactions figure (3a) 
197 Au(a!, rt) 200 Tl, figure (3b) 197 Au(a, 2n) 199 Tl, figure (3c) 197 Au(a,3n) 198w '*Tl, 
figure (3d) 197 Au(a,4n) 197 Tl and figure (3e) 197 Au(o;,/?3n) 197m Hg are summarized and 
plotted as a function of a-particle bombarding energy. The absolute errors which is 
generally = 8% consists of uncertainties due to target foil thickness (2%), the beam 



(a) 



197-Au(a. n)200-Tl 




30 40 50 60 

ENERGY IN MeV 



(b) 



197-Au(a, 2n)199-Tl 
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Figure 3(a-d). 



(d) 197-Au(a,4n)197-Tl 




40 50 60 

ENERGY IN MeV 



/ n i en ivr 



(e) 



197-Au(a, p3n)197m-Hg 




30 40 50 60 

ENERGY IN MEV 

Figure 3a-e. The total residual production cross-section in mb for the reaction (a) 
197 Au(a, n) 200 Tl (.), (b) 197 Au(a, 2n) I99 Tl (.), (c) 197 Au(a, 3n) 198m ' fi Tl (= in- 
state, a =* ^-state), (d) 197 Au(ct, 4n) 197 Tl (a) and (e) 197 Au(a, p3n) 197w Hg () for in- 
state) are plotted as a function of a-particle bombarding energy. Solid lines are the 
geometry dependent hybrid model and the -- or -o- lines are the hybrid model fits 
using the code Alice/91 with N = 4, N n = 2,N P = 2 and N h = 0. 



current integration (2%), the detector efficiency (5%) and the analysis of the 7-ray 
spectra (statistical uncertainty) generally (< 2%). The uncertainties caused by the large 
size of the irradiation area and the non-uniformities of the target contribute about 5% to 
the average error of the cross-section. The uncertainty given for the energy values are 
those of target thickness only and beam energy resolution was = 200.0 keV only. 

4. Comparison with theoretical predictions 

We have investigated the nuclear mechanism of the nuclear reaction induced by a- 
particles using the computer codes Alice/91 [17] and Stapre [18]. The Alice/91 code 
[17] describes the process of equilibrium evaporation of particles and 7-rays in terms of 
the Weisskopf and Ewing model [24] and pre-equilibrium reaction mechanism according 
to the hybrid and geometry dependent hybrid model [1,2]. Whereas in the Stapre code 
[18], the evaporation of particles and 7-rays are treated in the framework of the statistical 
model with consideration of angular momentum and parity using Hauser-Feshbach 
formalism [25]. For emission of the first particle pre-equilibrium decay is also taken into 
account. The pre-equilibrium emission of the particles was treated in the Stapre code [18] 
in the framework of exciton model [3]. The high energy parts of the excitation functions 
are dominated by pre-equilibrium reaction mechanism whereas the low energy parts are 
dominated by evaporation with its characteristic peak. 

Most of the nuclear reaction models [1-6, 26, 27] to treat the pre-equilibrium phase of 
reactions leading to the formation of a compound nucleus are semi-classical in nature and 
have been used with considerable success in describing experimental data pertaining to 
the equilibration process, mainly the forward peaked hard component observed in the 
continuous spectra of light ejectiles and the high energy tails seen in the excitation 
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Excitation function and isomeric cross-section ratios 

functions of activation cross-sections. However, most of these models employ one or both 
of the two basic concepts, the intranuclear cascade model [INC] [4,5] and Griffin's 
statistical model of intermediate structure [3] (SMIS). The integral excitation functions of 
a-induced reactions have been discussed by several authors [7-10] considering models of 
the compound nucleus as well as pre-equilibrium reactions, as mentioned above. They 
conclude that the theory of pre-equilibrium reactions is helpful in explaining the 
mechanism of a-induced reactions. The hybrid model for pre-equilibrium reaction was 
proposed by Blann [1] which provides in some way a marriage between the simple SMIS 
model of Griffin [3] and the INC model using the more elaborate master equation 
approach due to Harp et al [4,5]. However, there were some shortcomings in the hybrid 
model. It gives incorrect cross-sections and spectral distributions when applied to nucleon 
induced reactions at medium energies. Comparison with the INC model indicated that in 
the hybrid model, deficiencies resulted from the failure to take properly the enhanced 
emission from the nuclear surface. These deficiencies were partly rectified by 
reformulating the hybrid model as a sum of contributions, one term for each entrance 
channel impact parameter. In this way the diffuse surface properties sampled by the 
higher impact parameters were crudely incorporated into the precompound decay 
formalism in geometry dependent hybrid model (GDH) by Blann [2]. In the present work 
our excitation functions are calculated on the basis of hybrid model [1] and geometry 
dependent hybrid model [2] using the program Alice/91 [17] on the ND-500 and Super- 
32 computer at our centre. 

The statistical part of the Alice/91 [17] can account for a large variety of reaction 
types. The evaporation of neutrons, protons and clusters such as deuteron and a-particles 
was considered [24]. The binding energies and <2- values use d in the present code were 
based on experimental masses. The Alice/91 code [17] stores experimental masses in a 
'data file. Whenever the nuclear masses are not tabulated, these were calculated from the 
Myers and Swiatecki mass formula [28], liquid drop masses with pairing. We have used 
the option in the default version of GDH whereby only the first collision is localized 
according to the impact parameter [29] with all the higher order precompound terms 
being treated by the hybrid model i.e. using nuclear densities averaged over the nucleus 
and independent of the impact parameter. This is reasonable because the excitons can 
sample nearly the entire nuclear volume after a single scattering since the mean free path 
(mfp) values are = 4.74 fm. The inverse cross- sections were calculated using the optical 
model subroutine of Alice/91 [17], where the optical model parameters were those of 
Becchetti and Greenlees [30]. The Fermi level density used is of the form 



p(u) = (v^/12) (u - 

where u is residual nucleus excitation, a is the level density parameter taken as 
A/9MeV -1 which is the default option of the code and 6= ll/\/AMeV the pairing 
energy shift, with either a backshifted or standard pairing shift option. We have used the 
standard option. 

In the a priori formulation of the hybrid and geometry dependent hybrid model, the 
intra-nuclear transition rates are calculated either from the imaginary part of the optical 
model or from the free nucleon-nucleon scattering cross-section [31]. For particle 
energies exceeding 55 MeV the optical model parameters of Becchetti and Greenlees [30] 
used in the code are not applicable and thus at higher energies the calculation of the mean 



free path for intra-nuclear transitions should be calculated from nucleon-nucleon 
scattering cross-sections. Since our bombarding energy range was below 55MeV, we 
have used only optical model option for calculating the intra-nuclear transition rates. 

The Stapre code [18] has been used to calculate the theoretical cross-sections and 
isomeric cross- section ratios for a- induced reactions taking into account up to six 
sequential evaporation of particles and 7-rays. Each evaporation step is treated within the 
framework of the statistical model with consideration of angular momentum and parity 
using Hauser-Feshbach formalism [25]. For the emission of the first particle pre- 
equilibrium decay is also taken into account. Direct interactions were not considered. 

The parameters used in the code Stapre [18] were generally accepted ones. For 
calculation of the transmission coefficients of various particles for the code Stapre [18], 
the default global set of optical model parameters of Alice/91 code [17] was used. The 
transmission coefficients for 7-rays with transition energy 7 are expressed by the 7-ray 
strength function /XL (7) for the multipole radiation of type XL. For the El strength 
function, the Brink-Axel [32] model with global parameters was used and for Ml, E2, 
M2, E3 and M3 radiations, the Weisskopf model [33] was used. For energies, spins and 
parities of the discrete levels of the residual nuclei, the lowest 15-20 levels in [22] were 
used. The level density formalism of the back-shifted Fermi gas model expressed by Lang 
[34] was used for the continuum excitation energy region. The level density parameter 
a = A/9 was used for all nuclei, where A is the mass number of the nucleus. The spin 
distribution of the level density was characterized by the effective moment of inertia O e ff 
or better by its ratio to rigid body moment of inertia 6 r j g (77 = 6 e ff/6 r i g ). Since isomeric 
cross-section ratios are expected to depend strongly on the effective moment of inertia, all 
the calculations were performed for 77 = 0.5 and 77 = 1.0. 

In the Stapre code the pre-equilibrium emission of the particles were treated in the 
framework of exciton model [3] having the following ingredients. For the initial exciton 
configurations (PQ/ZO) we used (4,0) for a-particles [17]. The transition rates were 
calculated by the formulae of Williams and Cline model [35]. The average residual two- 
body matrix element that appears in the transition rates A+,A and A_ formulae as 
proposed by Kalbach-Cline [36] is expressed as 

|M| 2 = FM A~ 3 E~\ (5) 

where E is the excitation energy of the composite system. The quantity FM is a constant 
with the dimension of (MeV) 3 and generally treated as free parameter so as to get a good 
fit to the experimental data. The values used for the a-particle induced reactions were in 
the range FM = 750- 1000 MeV 3 . 

In Alice [17] code, we have taken an initial exciton configuration n = 4 
(n n = 2, n p = 2, rih = 0) which is equivalent to a break-up of the incoming a-particle 
in the field of the nucleus and the nucleons occupying excited states above the Fermi 
energy gives a better description of the excitation function compared to other configura- 
tions for the a-particle bombarding energies up to 55.0 MeV. Therefore we have mostly 
used no = 4 (n n = 2, n p = 2, n/, = 0) configuration in most of our calculations. The total 
production cross-section in mb for the reactions in figure (la) 193 Ir(a, 2n) 195 Au (), 
figure (Ib) 193 Ir(a, 3n) 194 Au (), figure (Ic) ( 193 )> 191 Ir[a, (4n)2n] I93 Au [(), D ] and 
figure (Id) 191 Ir(a, 3n) 192 Au () are plotted as a function of a-particle bombarding 
energy. The solid line is the geometry dependent hybrid model and the circle-line (-o-) is 
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the hybrid model fits using the code Alice/91 with NO = 4,N n = 2,N P = 2 and N h = 0. 
In figure 2a the total residual production cross-section in mb for the reaction 
181 Ta(o;, n) 184OTlS Re [symbols => g, n => m] are plotted as a function of a-particle 
bombarding energy. The solid line is the geometry dependent hybrid model fit using the 
code Alice/91 with N = 4,W n = 2,N P = 2 and N h = 0. The -*- and -o- lines are 
Stapre code calculations with FM = 1000 (MeV) 3 and FM = 750 (MeV) 3 respectively. 
In figure 2b and table 3 the isomeric cross-sections and their ratios are presented. The 
isomeric cross-section is practically constant at low bombarding energy with some 
tendency to increase at higher bombarding energy which means that at higher bombard- 
ing energy there is some tendency to populate the isomeric state preferentially. However, 
at very high excitation energy both the states are populated equally as seen in table 4 for 
the reactions !85 Re(a, an) 1 4m ' s Re. In figure 2c the total residual production cross- 
sections in mb for the reaction 18l Ta(o;, 2n) 183 Re [ symbols] are plotted as a function of 
a-particle bombarding energy. The solid line is 'the geometry dependent hybrid model fit 
using the code Alice/91 with NO = 4,7V,, = 2, N p = 2 and Nf, = 0. In figures 3a~e the 
total residual production cross-section in mb for the reactions figure (3a) 197 Au(a, n) 200 Tl 
(), figure (3b) 197 Au(a, 2n) 199 Tl (.), figure (3c) 197 Au(a, 3n) 198 *' M T1 (B => m-state and 
D =g-state), figure (3d) 197 Au(a, 4n) 197 Tl () and figure (3e) I97 Au(a, p3n) 197m Hg 
(a = m-state) are plotted as a function of a-particle bombarding energy. The solid lines 
are the geometry dependent hybrid model and the circle-dashed lines (-o-) are hybrid 
model fits using the code Alice/91 with NO = 4, W,, = 2,N P = 2 and Nh = 0. 

In the energy range of the measurement both geometry dependent hybrid model and 
hybrid model fit all the excitation functions reasonably well taking limitations of the 
calculations into account. Considering the multitudes of uncertainties in pre-equilibrium 
calculations such as (i) range of equilibrium and pre-equilibrium reaction cross-section 
involved and (ii) in parameters such as inverse reaction cross-sections and level densities 
etc. Blann [31] considered that a result which is within a factor of two (i.e. 0.5 to 2.0) of 
the experimental result in absolute cross-section and which generally has the correct 
spectral shape and variation of yield with excitation energy is an encouraging result. 



5. Conclusion 

A consistent set of 1 8 excitation functions has been measured for a-induced reactions on 
Ir, Au, Ta and Re targerts. The reliability of the cross-section was checked by comparing 
our measurements on Au with previous measurements. As in the previous studies, we 
have also observed that the excitation functions of the reactions induced by a-particles 
show structure with an initial peak with a slowly decreasing high energy edge becoming a 
flat plateau. The lower energy peak corresponds to particle emission mainly by 
evaporation by completely equilibrated compound nucleus. Whereas the low cross- 
section higher energy flat plateau corresponds to pre-equilibrium emission of nucleons at 
initial stages of the nuclear reaction. Comparison of the experimental data with hybrid 
and GDH model calculation shows reasonably good agreement with GDH model fitting 
high energy shoulder better indicating geometry dependent (surface) effects are important 
at higher energies only. Therefore, we can conclude that the dominant mechanism for 
reaction induced by a-particle in the lower energy range of these measurements is fusion- 
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evaporation and only at the higher energy the pre-equilibrium emission becomes 
dominant. 
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Abstract. A self consistent formulation of the Jeans instability of a dusty plasma with proper 
inclusion of charge dynamics is described. It is shown that charge fluctuations significantly affect 
the Jeans as well as the Buneman mode. For plasma particles (electrons and ions) in local thermal 
equilibrium, the Jeans length Aj is given by Aj X S F(R, e, /3/rj), where A g is the Debye length of 
the charged grains, R is the square of the ratio of the Jeans to the plasma frequency of the grains, is 
the square of the ratio of the Debye length of the grains and the plasma particles and (3/rj is the ratio 
of the attachment to the decay frequency of the electronic charges to the grain surface. The 
functional form of F is given in the text. Numerical investigation of the Jeans-Buneman mode for a 
two and three component plasma shows that the Jeans mode dominates at Ao <C 1 (where k is the 
wave number and AD is the Debye length of plasma particles), whereas at kX ^> 1 only the 
Bunernan mode operates. Charge fluctuations reduce the area of overlap of the two modes. 
Furthermore, in the absence of gravity, there exists a new, charge fluctuation induced unstable mode 
in a streaming dusty plasma. Astrophysical applications of the results are discussed. 

Keywords. Jeans-Buneman instability; Jeans length; charge fluctuations. 
PACSNos 52.25; 52.35 



1. Introduction 

Dust grains are present everywhere in the universe. Dark haloes in galaxies [1], infrared 
emission from circumstellar shells of early as well as late type stars [2,3], spectral 
features of the emission from active galactic nuclei (AGN) [4], ring and spoke formation 
in planetary systems [5], structure of cometary tails [6], are all attributed to the presence 
of micron sized grains. Back on earth, grains are present in uncontrolled (e.g. in the 
earth's magnetosphere and ionosphere) as well as controlled (e.g. Tokamak plasmas, 
plasma arcs etc.) environments. The grains are generally composed of graphite, silicate 
and metallic compounds in cosmic environments. The size distribution of the grain has 
been investigated by comparing the observed interstellar extinction curve with the 
theoretical one. The observed extinction curve do not fit to a single size and single 
composition implying that the sizes, masses and electric charges of the grains vary in a 
wide range. For example, grain sizes may vary between macroscopic (few cm) and 
microscopic (10~ 8 cm) scales, masses vary between 10~ 5 -10~ 15 g and electronic charges 
between (1 10 6 ). Charging of a grain takes place because of the radiative and plasma 
background in which it is immersed in most of the astrophysical situations. It is quite 
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possible that the sizes and masses of the grains might depend on the same plasma 
environment, e.g. Coulombic coagulation [7] could lead to bigger and heavier grains 
while breaking of the massive grains due to Maxwellian stress could be the source of 
lighter and smaller grains. These charged grains couple with the plasma dynamics 
through electromagnetic (EM) forces. A dusty /grainy plasma is thus a three component 
plasma consisting of electrons, ions and charged grains. 

Collective behaviour of a dusty plasma has been investigated in considerable detail in 
the recent past [5,8-11]. As most of these studies pertain to fixed charge on the grain, 
waves and instabilities in such a plasma are not very different from that of three 
component ion-electron plasmas. It is through the inclusion of charge dynamics [12-19] 
or mass dynamics [20] that new collective features have emerged. 

The gravito-electromagnetic coupling of an atomic plasma in local thermal equi- 
librium was studied in the context of star formation by Eddington [21]. He showed 
that self-gravity is unimportant for such a plasma as R = (Jj/u^ g ~10~ 36 , where wj = 
(47rGn g m g ) 1 ' 2 is the Jeans frequency and cj pg = (47rn g <2 2 //n g ) is the plasma frequency 
of the grain with G as the gravitational constant and n g , m g and Q as the number density, 
mass and charge of the grain respectively. Therefore, the scales at which the two forces 
operate are widely separated. The formation of large scale structures have entirely been 
attributed to gravitational condensation of matter whereas the radiation processes are 
attributed to electromagnetic interactions of the plasma particles. However, for typically 
cm sized grains (m g ~ 10~ 5 gm), 7?~0(1) and thus the two scales may overlap. The broad 
spectrum of observed masses and charges [22, 23] suggests that the dynamics of a dusty 
plasma can be studied in any of the following regimes: (a) EM GF (gravitational 
force), (b) EM force ~ GF and (c) EM force <C GF. The first case corresponds to the 
plasma processes like radiation, heating etc., the second case is thought to be the cause of 
Spoke formation in Saturn's ring, thickness of the Jovian rings etc., and the last case 
corresponds to the formation of stars, clusters etc. [24]. 

An important distinction should be pointed out between the R ~ 1 and Fgrav ~ ^Em 
(Fgrav and F Em are the gravitational and the electromagnetic force on the grains). Whereas 
R ~ 1 leads to unusual collective behaviour of a self-gravitating dusty plasma due to 
overlap of two scales, Fgr av ~ FEW implies merely the balance of the two forces whose 
origin could be quite different. For example in the planetary rings, the self repulsion of 
the grain is balanced by the gravitational attraction of the planets for micron and 
submicron sized grains [25]. 

Recently, the Jeans instability of a dusty plasma has been studied by several authors 
[26-31]. Time dependent nonlinear solutions of a self-gravitating three component dusty 
plasma showed [27] that when self-gravity is annuled by the electric polarization of a 
collapsing plasma i.e., R ~ 1, plasma may become gravitationally unstable due to the 
properties of the background plasma. Therefore, merely balancing the self attraction of 
the grains with self repulsion is not sufficient for determining the fate of condensation. 
Avinash and Shukla [29] studied the Jeans instability of a dusty plasma with the inclusion 
of ion dynamics and found enhanced gravitational condensation of the grains. Later, 
Pandey and Dwivedi [30] showed that the proper inclusion of ion dynamics does not 
destabilize the Jeans mode any further but the collapse of the grain follows the same 
dynamics as the collapse of the neutral matter under gravity. This is attributed to the fact 
that ions are responding at the free fall time scale of the grains and as a result they shield 
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/vii me previous siuuies, wim me exception or me last one, exclude cnarge dynamics 
from consideration. However, the grain charge is a dynamical variable and in the present 
work we investigate the Jeans instability with proper charge dynamics of the grains. We 
consider an infinite, homogeneous, unmagnetized, nonrotating plasma and, invoking 
'Jeans swindle', ignore the zeroth order gravitational field. However, as has been 
demonstrated [27], the construction of a proper equilibrium does not change the plasma 
dispersion relation, as one can always consider homogeneity in a limiting sense. We first 
consider the dusty plasma with electrons and ions in local thermal equilibrium and find 
that in the presence of charge fluctuations, the Jeans length gets modified, implying that 
the charge dynamics has a bearing on the gravitational collapse of the matter. 

Next, we consider the gravitational instability of a two component dusty plasma with 
streaming ions. We find that when charge fluctuations are absent, a new gravito- 
streaming mode may give rise to the accretion of matter, when self gravity is annulled by 
the thermal pressure and self repulsion of the grain. Inclusion of charge fluctuations 
modifies both gravitational as well as streaming modes. Also, a new charge fluctuation 
induced instability in the presence of streaming ions exists in such a plasma. However, the 
growth rate of this new mode is much smaller than that of the Buneman mode. 

Lastly, the dynamics of a three component dusty plasma is analysed. In the presence of 
streaming ions, the picture is similar to the previous case i.e., the Jeans and the Buneman 
modes dominate in two different windows with a small overlap of the two. The streaming 
instability sets in at a smaller wavelength in the presence of charge fluctuation than in the 
absence of it. The growth rate of the Jeans mode is enhanced in the presence of charge 
dynamics of the grains. Also, a charge fluctuation induced instability exists in such a 
plasma with a growth rate much smaller than that of the Buneman mode. 

The paper is organized in the following fashion. In 2, basic equations are given and 
the equilibrium is discussed. Section 3 describes the linearized equations and gives the 
general dispersion relation. Also, a special case of the dispersion relation is discussed 
when electrons and ions are in local thermal equilibrium. Section 4 describes the Jeans- 
Buneman mode for two and three component dusty plasma systems. In 5 a summary of 
the results is given and possible applications to interstellar media are indicated. 



2. Basic equation 

A three component dusty plasma consisting of electrons, ions and charged dust grains, is 
assumed to have grains of equal radius and carrying identical charge due to the collisional 
attachment of the electrons and ions to the grain surface. Effects of the photoemission, 
secondary emission etc. is ignored in the present study. It should be mentioned here 
that this assumption is not entirely satisfactory for the thermal grains. The charging 
and the discharging of the grain give rise to sink and source terms in the continuity 
and momentum exchange equations. Owing to the large mass difference between the 
plasma and the grains, we assume that the gravitational potential is determined mainly 
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equations: 
Continuity equations: 

^+V-(n Q v Q ) = -/3 a (n tt -rt QO ), (1) 

r^ 

^+V-(n g v g ) = 0, (2) 

where a e,i\g = grain. The subscript means equilibrium quantities and (3 a = I a n & / 
q a n a is the attachment frequency of the plasma particles to the grains. 

Equations of motion: 

dv Q 



dv 



(2 g V0 m g ng\/i/}. (4) 

U.fr 

Charge equation: 

. dQ 
~ = / e (Q,0)+ 7/02,0), (5) 

where the currents are given by [32, 33] 

e exp - i , (6) 

L 1 e J 

-^] (7) 

and the two Poisson's equation are 

V 2 = 47re[n e - n,-] - 47rgn g , (8) 

V 2/ = 47rGrc g m g , (9) 




where d/dr = 3/5r + v- V and m, , u, 7, 0, ^ denote the mass, number density, velocity, 
temperature, electrostatic and gravitational potential respectively. The quantity g - is 
the potential difference between the grain surface potential and the background plasma 
potential so that the unperturbed charge > is given by Q Q = C(0 g - 0) where C is the 
capacitance of the grain and is given by [34], C = aexp (-a/A D ) with a as the radius of 
the grain. The exponential factor reflects the screening of the grain charge by the plasma 
particles. In equilibrium, electron and ion currents equal each other and thus the grain is 
at the 'floating potential' . Note that for n o/;o ~ 1 (a situation relevant to most of the 
planetary system [10] and interstellar media [9]), /3 e = ft i.e., electrons and ions get 
attached to the grain at the same rate after the initial build up of a negative charge on the 
grain similar to the behaviour of plasma particles near the sheath [35], The above set of 
equations (l)-(9) form the basic set of equations for the linear analysis. 
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guaranteed by the existence of opposite charges, there is no such thing as opposite mass 
and thus, there is no way to make the gravitational potential disappear in zeroth order. 
Therefore, the problem of a self gravitating system is in general an eigenvalue problem. 
However, important physical insight can be gained by invoking 'Jeans Swindle' i.e 
ignoring the zeroth order potential field of gravity. To reiterate, we consider the linear 
gravitational instability of an infinite, homogeneous plasma characterized by 

n e o = n i0 + Zrt g o, -00 = </>o = 0, v/ = VQX, v eG = v g o = 0, 
where VQ is the ion streaming velocity and Z = Q/e. The linearized set of equations are 

"' 



dt 



dt 



+ V-(n tt iV + n Q 
+ n g0 V-v g i =0, 



(10) 
(11) 



where the subscript 1 means perturbed quantities. 



n a o 



'gi _ 



dt 



dt 



m 



(12) 
(13) 

(14) 



where 77= (e\I e o\/C)(l/T e + I/WQ), wo = 7/ - e0/o with TI as the ion temperature and 0/ 
as the equilibrium floating potential of the grains and C a = (r tt /m Q ) ' is the thermal velo- 
city of the aih species and p a \ is obtained from the definition of /3 a by perturbing it around 
/3 Q o- Taking the spatial and temporal dependence of all perturbed quantities as 
exp [-i(ujt kx}], we get the following relation between fluctuation density and potentials: 



v 



+ 



-&& 



(15) 



(16) 



and 



m, 



fQ 

\m Zn g 
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where use has been made of Ai = A(( g i/ g o) x(<2i/<2)), Pel = A?((gi/go)+ 
x(<2i/<2)) with x = l/((?V<%o) - 1) [14]. Then plugging equations (15)-(17) in the 
Poisson's equation leads to the following dispersion relation: 



i ^ i .\ i \ . 

1 i 1 -I -- _ A + 1 H 



, 
tu +- ir u + w - 

where 



2 



(w + /A) 2 - k 2 Ct + ifrkVo + x/3f 

(20) 

where cD = w kVo is the Dopier shifted frequency, ( = 47rn a q a 2 /m a is the plasma 
frequency of the ath species and = n e o/io- Subscript has been dropped from the 
equilibrium quantities. The above dispersion relation in the limit uj pg > 0, k\ e <C 1 , 
reduces to equation (18) of Bhatt and Pandey [14]. First we analyse the gravitational 
instability of a dusty plasma when the plasma particles are in thermal equilibrium. For 
this case the thermal response of the plasma particles is much faster than the grain's 
response to the perturbations i.e., &C g <C w<C kCf < kC e where k is the wave number. 
Assuming j3 e = A = /? we get the following dispersion relation: 



2 



where A D = A e Aj/[A^ + A?] 1//2 is the effective plasma Debye length with X e ,i = 
(T et i/47rn et ie 2 ) ' . The dispersion relation (21) coincides with that of Mohanta et al 
[31] in the absence of an external magnetic field for a cold dust grain. Before undertaking 
a complete numerical analysis of the above dispersion relation, the effect of the charge 
fluctuation can be studied by solving equation (21) perturbatively by treating 77 ~ (3 < uj. 
Then the interesting root is 



(22) 

where Ag = T g /4-rm g Q^ is the Debye length of the grains. From the above expression it is 
clear that the attachment of the grains causes enhancement of the gravitational 
condensation for some < k < kj, where k] is the Jeans wave number at which the right 
hand side of equation (21) goes to zero in the absence of charge fluctuations, i.e., = 0. 
As the wave number crosses the Jeans wave number, the first term in equation (22) 
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Figure 1. Curve 1 corresponds to the perturbative solution for (3/r] ~ 1 and 
x2 = /3/u pg =0.1. Curves 2, 3, 4 and 5 shows the Jeans mode for ^1 = f]/^ PS = 0, 
x2 = 0; xl = 0.1, z2 = 0.1 and xl = 0.1, x2 = 1 and xl = 1, x2 = 0.1 respectively. 
Heree= 10" 6 , R = Q.l. 



becomes real and thus it describes the usual dust acoustic wave [36]. Around k = k] the 
perturbative analysis breaks down (see curve 1 of figure 1 ; the parameters for which the 
plot is drawn are given below). Qualitatively one can explain these features as follows: 
Charge fluctuations damp the dust acoustic waves. The rate of damping is maximum at 
certain values of k and reduces to zero for both higher and lower k values. This will cause 
the shift in the balance of force in favour of gravity for some range of k. As a result, 
gravitational instability becomes robust and the Jeans length reduces compared to the 
fixed charge case. 

Next, we analyse equation (21) numerically. We solve equation (21) for e Ag/A^ = 
lQ~ 6 ,r)/w p g = (Q 0.1) and (3/uj p& ~(Q l),R = Q.l, numbers representative of the 
interstellar medium [10]. From figure 1 (curves 2 and 3) we see that in the absence of 
charge fluctuation, the growth rate goes to zero at smaller k\ D than in the presence of 
charge fluctuation, implying that the condensation of the grain is facilitated by charge 
fluctuations. An increase in the attachment frequency (/3 ^> 77) leads to the enhancement 
of the growth rate (figure 1 , curve 4), since attachment of the plasma particles favourably 
reduces the repulsive potential well of the grain, facilitating the condensation. Similarly, 
an increase in 77 implies that more plasma particles are seen in the vicinity of a collapsing 
grain, repelling the incoming grain in the process (curve 5). Thus, the growth rate 
decreases. With the increase of R the growth rate will increase, as it is the measure of 
relative strength of the attraction to the repulsion. 

Finally we calculate the Jeans length from equation (21). Jeans length is the critical 
length at which the instability vanishes. Then solving equation (21) for the marginal 
stability condition, we get 



A == 



/w 



(23) 
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and 



1/2 



1/2 



where e = A/A^. For R ~ 1, 




e< 1. 



(24) 



(25) 



From the above limiting expression for F it is clear that when e ~ 1, grain charge 
fluctuation plays no role in the gravitational collapse and Aj = A g . Therefore, the 
gravitational condensation of the grain can be halted only by the repulsive electrostatic 
field for a fixed charge. However when e <C I Jeans length is considerably modified by 
the charge dynamics of the grain. Depending upon the value of /3/r] the Jeans length can 
decrease or increase. For e = 0.1 we see from figure 2 that Aj/A g initially decreases (in 
comparison with the value of Aj/A g in the absence of charge fluctuations) and then after 
/3/T7 = 5.5 it shoots up. Therefore, for grainy plasma whose Debye length is much smaller 
than the Debye length of the plasma particle, charge fluctuations can facilitate or oppose 
the collapse of the grain by setting up the electrostatic field in phase or out of phase with 
the gravitational field. Above limiting cases can be applicable to different regions in 
interplanetary and interstellar media. For example, it is known that dark molecular clouds 
(DMC) contain grains of different charges depending upon the ambient plasma tempera- 
ture [37]. If the ambient temperature is between 0-1 OK then the grain may have 



0.8 



0.6 - 
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12 



Figure 2. Normalized Jeans length is plotted against j3/-q for e = 0. 1 and R = 1 . The 
horizontal curve is the value of normalized Jeans length in the absence of charge 
fluctuation. 
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Jeans-Buneman instability 

1, 1 charge whereas for temperature around 100K a grain may have 10 electronic 
ges. As e may become 0(1) in DMC [9] the first of the two cases may be relevent. 
above analysis indicates that for grains charge flip-flopping between -1 , 0, 4- 1, i.e., 
77 or for grains containing 10 or more charges /3 > 77, the charge fluctuation is of no 
sequence to the gravitational condensation of the grains. In a planetary environment 
: 1 and typically Q ~ (10 3 - 10 6 )<? [5, 10]. Then the charge fluctuations can cause 
rease as well as increase in the Jeans length. As a result charge fluctuations can 
litate or inhibit the condensation of the grains depending on the value of (3/r]. 
refore we see that for A g -C AD gravitational condensation of the grain can be facili- 
d or halted by an electrostatic field much beyond the length suggested by the analysis 
'andey et al [27]; Pandey and Dwivedi [30] for a fixed charge. The preliminary nature 
he above results must be noted as they have been carried out for very idealized 
ations. For a final application of these results to relevant physical problems (viz. star 
nation, formation of the planetary system etc.), a model calculation must be done. 



feans-Buneman mode 

t, we shall consider a two component dusty plasma with streaming ions. This is an 
il situation where all the electrons are attached to the grains. This type of plasma may 
it in some of the regions of dense interstellar clouds, where Zn g may become 
iparable to n ( [38]. Dynamics of such a plasma has been considered [39] in the context 
L tearing instability in the Saturnian rings. Then the dispersion relation equation (18) 
ices to 



1 



UJ 2 




p& 


D 


uj 2 + uj 2 -k 2 C 2 


r ;/3 1 


22 


1+ ft 


JL J PS 


[ UJ + ir)\ 





(26) 



D = [(u) 2 + 2i#<D) - fc 2 C 2 - iA- 



ore numerically solving the above 5th order polynomial, we analyse equation (26) in 
absence of charge fluctuations. Then the dispersion relation becomes 



[u, 2 + u, 2 - k*Cl - o, 2 g ][(u> 



u; 2 ,)] = 



(27) 



in the absence of gravity we have the usual Buneman mode [40]. Here we look for a 
ie in the vicinity of kVo w (k 2 Cf -fo;^) 1 ^ 2 in the presence of gravity. The unstable 
ns-Buneman root is 



(28) 



vided ujj a; 2 . The significance of this mode underlies the fact that even when self 
action is completely annulled by self repulsion, the Jeans-Buneman mode may trigger 




1.5 - 




2.5 



3.5 



Figure 3. Curves 1, 2 and 3 are the Jeans mode with xl = Q,x2 = 0; x\ = 0, x2 = 1 
and xl = 1, xl = 1 respectively for ^ = 0.1. Curves 4, 5 and 6 shows Buneman mode 
for a = VQ/CI = 1.1, e = 1CT 4 , R = with xl = 0, x2 = 0; xl = 0, *2 = 1 and 
xl = 1, x2 = 1 respectively and curve 7 shows the Jeans-Buneman mode for xl = I, 
x2= 1 with/? = 0.1, a = 1.1. 



off a large scale condensation of the grains. The growth rate is given by 



T ~ 2.10 



c 



m a n a n ; 



1/2' 



m 



1/2 



1/3 



Hz, 



(29) 



where, n, m are measured in cgs units. The importance of the above mode can be gauged 
from the following example. It is known that the self-gravitating clouds of masses 
10 2 10 6 M Q , where M is the solar mass, contain many dark cores with number densi- 
ties in excess of 10 4 cm~ 3 and masses < 10 2 Mo [9]. In these dark cores, nearly one per 
cent of the mass is in dust form. In such cores, typically n g w 4 x 10~ 12 (<2/10~ 5 )nH 
which for hydrogen density tin ranging between (10 3 cm~ 3 - 10 7 cm~ 3 ) and cm sized 
grain becomes (lcm~ 3 10~ 4 cm~ 3 ). Calculating ion density from n,-/nH~10~~ 8 
(/7 H /10 5 cm" 3 )" 1 ^ , we get n,- ~ (1CT 2 cm" 3 10~ 4 ) cm"" 3 . Assuming ions to be ionized 
hydrogen we get F~10~ 3 . Thus we see that the Jeans-Buneman mode may play an 
important role in the formation of large scale structures. 

Next, we solve the dispersion relation (26) numerically. We see from figure 3 that at 
large scale lengths only the Jeans mode operates, whereas at small scale lengths the 
Buneman mode is excited. There is an overlap of the two modes. The growth rate of the 
Jeans mode is smaller than that of the streaming mode. Charge fluctuations affect both 
gravitational and Buneman modes. In the presence of charge fluctuation, the growth rate 
of the Jeans mode increases and drops to zero at large A;A D . The Buneman mode, on the 
other hand, starts operating at much shorter wavelengths in the presence of charge 
fluctuations. The area of overlap of the Jeans and Buneman modes shrinks in the presence 
of charge dynamics. At kX D <C 1 the gravitational instability dominates over the 
streaming instability. In the other limit, i.e., fcAo ^ 1, it is the streaming instability which 
destabilizes the grain plasma system. In the very narrow window Ao ~ 0(1], the two 
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0.5 1 1.5 



Figure 4. Curves 1, 2 show die charge fluctuation induced unstable mode for R = 0, 
a = 1.1, e = 10~ 6 with jcl = 0,x2 = 1 and xl = I,jc2 = 1 respectively. 




15 20 

k\ D 

Figure 5. Curves 1, 2 and 3 shows Jeans mode for a three component plasma with 
xl = 0, x2 = 0; xl = 0, xl = 1 and xl = 1, xl - 1 respectively for R = 0.1. Here 
e = if)' 6 , 6 = 0.01, m/Zn s = 0.99, \J/X 2 e = 10~ 4 , wj,/u g = 10 3 , wj,/uj g = 10 2 . 



les overlap. In the absence of gravity, in addition to the Buneman mode, we find a 
charge fluctuation induced unstable mode existing in a streaming plasma as shown in 
re 4. The growth rate of this mode is smaller than that of the Buneman mode. This 
le appears analogous to the negative energy mode reported earlier by Jana et al [13]. 
sically as a = VQ/CI > 1, the ion wave can be viewed as a negative energy wave 
n the ion frame, so that dissipative effects associated with nonstreaming grains can 
5 rise to its amplification. 

lext, we study the stability of a three component plasma. We solve the dispersion 
tion (18) numerically for Boltzmannian as well as non-Boltzmannian electrons. The 
usion of electron dynamics does not affect the result significantly and the change in 
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Figure 6. Curves 1 , 2 and 3 shows the Buneman mode for R = with x\ = 0, x2 = 0; 
xl =0, x2= 1 and xl = 1, x2= 1 respectively. Here e= 10~ 6 , 6= 0.01, n,/Zh g = 0.99, 
X = 10~ 6 , Af/A2 = 10- 4 , cj2 fi /o;J g =10 3 , cjj./c^g = 10 2 . Curve 4 depicts Jeans- 
Buneman mode for R = 0.1, a = 1.1. 



the result is only 1%. In the parameter regime studied by us, last term in the numerator of 
equation (18) is unimportant. However, there exists the parameter window in which 
contribution of this term is comparable with the first two terms of the numerator. 

In figure 5 we plot the Jeans mode with and without charge fluctuations. As is evident 
from the figure, the modification due to charge fluctuations in the Jeans mode for a three 
component plasma is similar to a two component plasma. Figure 6 depicts the Buneman 
mode and also the Jeans-Buneman mode. The qualitative features are similar to the 
previous mode. We find that charge fluctuation in the presence of streaming ions also gives 
rise to negative energy mode. The behaviour of this mode is similar to that described in 
figure 4 for a two component plasma. We mention here that this result is at variance with 
the result of Bhatt and Pandey [14] who found perturbatively that the parameter regime for 
negative energy mode shrinks to zero in a three component plasma when 8 w 1 . 

5. Conclusion 

We have analysed the Jeans instability of a dusty plasma self consistently by properly 
considering the charge dynamics. It has been shown that charge fluctuations cause a 
modification of the Jeans length. The charge fluctuation can facilitate or hinder the 
condensation of the grains depending upon the value of /3/r). Implication of this result on 
Spoke formation in the Saturn's ring or on the thickness of the Jovian rings may be 
profound. For example, as the charge fluctuations modifiy the Jeans mode, the lifetime of 
self-levitation may be finite before the grains finally succumb to their self-gravity. In the 
absence of charge fluctuations the levitated grains will forever be hanging above the ring 
plane of the planets. In the interstellar media, even highly charged grains (e.g. ~ lOOe in 
the HII region) will undergo an enhanced condensation due to the dynamic nature of the 
grain charge. Therefore, our analysis indicates that charge fluctuations might help the 
gravitational condensation of the charged grains in different cosmic environment. 
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Analysis of two and three component dusty plasma do not qualitatively show different 
aviour. The Buneman and the Jeans modes both act in quite distinct zone with a very 
ow overlap region in the presence of charge fluctuations. For a two component plasma, 
find an analytical expression for the Jeans-Bunernan mode. The growth rate of 
aming instability is always higher than the Jeans instability. Therefore, the existence of 
ixed Jeans-Buneman mode will cause faster condensation of the grainy plasma than 
purely gravitational instability. A new charge fluctuation induced unstable mode is 
id from the numerical solution of the dispersion relation. It might be possible that when 
Buneman mode saturates, the charge fluctuation induced mode will be of significance, 
lased on the above analysis one can draw a general picture of the collapse of plasma 
ter. It might be possible that the collapse of the plasma matter may undergo several 
;es: Far away from the gravitating centre, gravitational attraction may trigger off a 
apse. The collapse may continue up to a certain scale length (say) LI at which stage 
;ma particles (namely ions) may develop a considerable flow (a > 1) which results in 
enhancement of the collapse rate. In the final stage which will be reached after a 
:ain scale length L, L <C L\ streaming ions go on accreting to the grains solely due to 
r flow. Thus, very close to the gravitating centre, its a purely Buneman mode which 
litates the condensation of the matter. 

t is known that massive objects surrounded by a charged plasma cloud (e.g. HJJ 
ion) will tend to accrete matter. The resultant flow of the plasma can be of importance 
stars in close binary systems [32], for compact objects (white dwarfs, neutron stars, 
hole etc.) which have been proposed as the possible energy source for the radiation 
ling out from a very localized region of the sky [41]. The accretion rate is calculated 
ely by the gravitational field. However, as one sees from the above analysis, the 
retion rate will become a function of the scale length. Very close to the gravitating 
tre the Jeans mode might be superseded by a Buneman mode and one gets an 
anced growth rate. The presence of a stellar wind in the interstellar medium supports 
h a picture. However, our analysis is for a very simplified model and thus the results 
suggestive rather than compelling. 

^.s a variety of plasmas is found in interstellar media, we believe that our results have 
ring on star formation, cluster formation etc. Charge fluctuations affect the 
idensation of the grains and thus, it might be possible that the rate of star formation 
the galaxies are different than the presently predicted rates. Further work in this 
;ction is required. 
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